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Abstract 



^ Wc study minimal implementations of Majorana and Dirac neutrino scenarios in F- 

^ theory GUT models. In both cases the mass scale of the neutrinos rrii, ~ M^^^^/ A^y arises 
^ from integrating out Kaluza-Klein modes, where Auv is close to the GUT scale. The partic- 
^ ipation of non-holomorphic Kaluza-Klein mode wave functions dilutes the mass hierarchy in 
comparison to the quark and charged lepton sectors, in agreement with experimentally mea- 
^ sured mass splittings. The neutrinos are predicted to exhibit a "normal" mass hierarchy, 
O with masses (ma, m2, mi) ~ .05 x (1, Oq^j,, aouT) eV. When the interactions of the neutrino 

> and charged lepton sectors geometrically unify, the neutrino mixing matrix exhibits a mild 
^ hierarchical structure such that the mixing angles 623 and 612 are large and comparable, 

^ while is expected to be smaller and close to the Cabibbo angle: ^13 ~ 6*^ ~ ckqut ~ O-^- 
This suggests that ^13 should be near the current experimental upper bound. 
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1 Introduction 

The observation of neutrino oscillations (T||^ has revealed that neutrinos have small non- 
zero masses. However, non-zero neutrino masses cannot be accommodated in the Standard 
Model without introducing extra ingredients. As such, neutrino physics offers a concrete 
and exciting window into physics beyond the Standard Model. 

The seesaw mechanism is perhaps the simplest theoretical model which describes small 
neutrino masses. By introducing very heavy right-handed Majorana neutrinos, the seesaw 
mechanism produces an effective light Majorana mass for the left-handed neutrinos. For 
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the masses of the left-handed neutrinos to be consistent with experimental bounds, the 
right-handed neutrinos must have Majorana masses around the scale Auv ~ 10^^ — 10^^ 
GeV, which is close to the GUT scale. Hence, the seesaw mechanism suggests that neutrino 
physics should be somehow related to the dynamics of GUT theories. 

However, in four-dimensional GUT models additional ingredients must be added just to 
accommodate the seesaw mechanism. For instance, in SO{10) GUTs, this necessitates ad- 
ditional fields transforming in higher dimensional representations developing suitably large 
vevs, or higher dimension operators (see for example |3] for a review of such mechanisms in 
the context of four-dimensional GUTs). Therefore, it is worth asking whether string theory 
may offer new insights into neutrino physics. 

In recent work on GUTs realized in F-theory (F-theory GUTs) the observation that 
Mgut/Mpi ~ 10~^ is a small number has been promoted in [4||5] to the vacuum selection 
criterion that there exists a limit in the compactification where it is in principle possible 
to decouple the effects of gravity by taking Mpi — )■ oo, with Mqut kept finite. See j6-21 



for some other recent work on F-theory GUTs. Aspects of flavor physics in F-theory GUTs 
have been studied in [12|, where it was shown that with the minimal number of geometric 
ingredients necessary for achieving one heavy generation, the resulting flavor hierarchies 
in the quark and charged lepton sectors are in accord with observation. The aim of this 
paper is to extend this minimal framework to include a neutrino sector with viable flavor 
physic^ 

We study both Majorana and Dirac neutrinos in minimal SU{^) F-theory GUTs, finding 
scenarios which lead to phenomenologically consistent models of neutrino flavor. In both 
cases, integrating out massive Kaluza-Klein modes generates higher dimension operators 
which lead to viable neutrino masses. The neutrino mass scale rriv is roughly related to the 
weak scale and a scale close to Mqut through the numerology of the seesaw mechanism: 

^ _ weak /-j^ -|^\ 



A 



uv 



In the Majorana scenario, an infinite tower of massive modes trapped on a Riemann surface 
play the role of right-handed neutrinos, and generate the F-term 

through an effective Kaluza-Klein seesaw mechanism. When Hu develops a vev {Hu) ~ 



^See 22 for other forthcoming work on flavor physics in the context of F-theory GUT models. 
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-^weak this induces a Majorana mass. In the Dirac scenario, the D-term 



Auv 



^Dirac I d^^^l ^ (1.3) 



is generated by integrating out massive modes on the Higgs curve. Supersymmetry breaking 
leads to an F-term for of order Fh^^ ~ fiHu ~ ^weak which induces a Dirac mass. We 
show that the participation of an infinite tower of massive states can boost the overall scale 



of the neutrino masses. This is welcome, since the two higher dimension operators (1.2) 



and (1.3) with scale Auv = Mqut would produce light neutrino masses which are slightly 
too low. 

Owing to the rigid structure present in F-theory GUTs, it is perhaps not surprising that 
the supersymmetry breaking sector of [o] naturally enters the discussion of neutrino physics. 
In [9] , the absence of a bare fx term in the low energy theory was ascribed to the presence of 



a U{1) Peccei-Quinn symmetry, derived from an underlying Eq GUT structure. This choice 



of U{l)pQ charges turns out to also exclude the higher dimension operator (1.2) appearing 
in the Majorana scenario. Interestingly, we find a unique alternative choice of U{1) charge 
assignments which is simultaneously compatible with a higher unification structure and the 
operator ([L2|)p1 



Estimating the form of the Yukawa matrices for the two operators (1.2) and (1.3), 
we find that in both scenarios the neutrinos exhibit a "normal" hierarchy, where the two 
lightest neutrinos are close in mass. The participation of Kaluza-Klein modes dilutes the 
mass hierarchy in comparison to the quark and charged lepton sectors. More precisely, the 
resulting neutrino mass hierarchy is roughly: 

1 /2 

mi : 7712 : mg ~ acuT ■ c^gut '■ 1 (1-4) 

which is in reasonable accord with the observed neutrino mass splittings. 

The structure of the neutrino mixing matrix depends on whether the neutrino and 
lepton interactions localize near each other, or are far apart. When these interactions 
are geometrically unified at a single point, the mixing matrix displays a mild hierarchical 
structure. The two mixing angles 6*12 and 623 are found to be comparable, and in rough 
agreement with experiments. The mixing angle ^13, which measures mixing between the 
heaviest and lightest neutrino (in our normal hierarchy), is predicted to be roughly given 
(in radians) by: 

^13 ~ ~ «GaT ~ 0.2, (1.5) 



^Even though this new U{1)pq does not change the general scenario of F-thcory GUTs, it does change 
some of the detailed numerical estimates of the "PQ deformation" away from minimal gauge mediation 



studied in [9 20 . It would be worth investigating this further. 
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where 9c denotes the Cabibbo angle. These results, in conjunction with the analysis of 12 



points towards the possibility of a higher unification structure. Along these lines, in both 
the Majorana and Dirac scenarios we present models where all of the interactions of the 
MSSM unify at a single interaction point in the geometry. 

We also study geometries where the neutrino and lepton interaction terms do not unify. 
In this case, the neutrino mixing matrix is a generic unitary matrix with no particular 
structure. As a result, large mixing angles are expected, and in particular the angle ^13 
should be close to the current experimental upper bound. Assuming that the neutrino 
mixing matrix is given by a random unitary matrix, we explain how randomness suggests 
that 9i2 and ^23 should be comparable, while ^13 should be slightly smaller, which is in 
qualitative agreements with neutrino oscillation experiments. 

The organization of the rest of the paper is as follows. In section |2] we review the main 
features of neutrino physics. Section [3] provides a short review of those aspects of F-theory 
GUTs which are of relevance to neutrino physics. We present a minimal implementation 
of the Majorana scenario in section |4j In this same section, we study the presence of 
monodromies in seven-brane configurations, and explain the crucial role this geometric 
ingredient plays in the Kaluza-Klein seesaw. In section[5]we estimate the Majorana scenario 
Yukawas. Next, in section [6] we discuss a minimal Dirac mass scenario, which surprisingly 
exhibits similar numerology to that of the Majorana scenario. Our results for the neutrino 
masses and mixing angles are compared with experiments in section [7} Section [8] contains 
our conclusions. Appendices A, B and C discusses other aspects of F-theory neutrinos, and 
Appendix D contains a discussion of probability measures for random unitary matrices. 



2 Review of Neutrino Physics 

In this section we review the main features of neutrino physics. We first describe background 



material on the masses and mixing angles of the neutrino sector in subsection 2.1, and then 



review current observational constraints in subsection 12.21 This is followed in subsection 



2.3 by a brief discussion of the suggestive appearance of UV physics in the neutrino sector 



and potential sources of tension with string based models. 
2.1 Neutrino Masses and Mixing Angles 

In this subsection we define the neutrino masses and mixing angles. In order to maintain 
continuity with the superfield notation employed later, we let L denote the lepton SU (2) 
doublet superfield of the MSSM, and the left-handed neutrino component of this doublet. 
We shall also denote by the charged lepton component of the doublet L, and by Er the 
right-handed charged lepton superfields. We emphasize that this notation is adopted for 
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notational expediency. Indeed, at the energy scales where the neutrinos develop masses, 
supersymmetry has already been broken. 

Neutrino mass can in principle originate from one of two possible effective chiral cou- 
plings, which below the electroweak symmetry breaking scale can be written as: 

WM.ior.n.^mp-NlNl (2.1) 
W^Dirac D mf^' ■ NlNl (2.2) 

where in the second case, denotes a right-handed neutrino, and i,j = e,fi,T index the 
three generations of left-handed neutrinos. These mass terms correspond respectively to 
Majorana and Dirac mass terms. The full lepton sector of the theory can then be written 
as: 

W^Lepton D ml;) ■ NIN^ + rr^ ■ E^E^ (2.3) 



where the first term corresponds to either of the two mass terms given in lines ( 2.1[ ) and 



(12J. 

As usual, we introduce matrices U^^ and and matrices [/|''* and diagonalizing 
the mass matrices in the lepton sector: 

Ui^'^m^'''^ (ul^^y = diag(mi,m2,m3) (2.4) 
f/f m(') (uj^y = diag(me, m^,m^). (2.5) 

Using the fact that N}^ and E}^ transform as SU (2) doublets, we can define a mixing matrix. 



as in the quark sector. The neutrino mixing matrix is given by 23 , 24 





Ue2 


Ue3 








Url 


Ur2 


Ur3 



UpMNs = ui^Uu^'A^ = \ \. (2.6) 



Introducing the parametrization of the unitary matrix in terms of the mixing angles 
< 9ij < 90°, we can write: 

(C12C13 ■5i2Ci3 SisC '"^ \ 

-S12C23 - Ci2S23-5l3e*'^ C12C23 - Si2S23Si3e*'' S23C13 • ^a, (2.7) 
iS iS I 

S12S23 " Ci2C23Si3e — C12S23 — Si2C23Si3e C23C13 / 

where Da = diag(e*"i/^, e*"^/^^ X)^ q.^ = cosOij and Sij = sinOij. Here, 5, ai and 0^2 are CP 
violating phases. In the Dirac scenario, only 6 corresponds to a physical phase, whereas in 
the Majorana scenario all three angles are physical. 
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2.2 Experimental Constraints 



Neutrino oscillation experiments have established that neutrinos are indeed massive [T| 
[2]. While we do not know the absolute mass eigenvalues mi, m2 and ms, experiments 
have measured small mass splittings. It is important to note that neither the relative 
spacing between the three neutrino masses, nor the lower bound on the neutrino masses has 
been established. There are three relative mass spacings which are in principle possible, 
corresponding to mi ~ m2 ~ ms, mi < m2 ^ ms and ms ^ mi < m2, which are 
respectively known as degenerate/democratic, normal hierarchy and inverted hierarchy mass 
spectra. As reviewed in (25|[26] , solar and atmospheric measurements of neutrino oscillation 
lead to the mass splittings: 

Am^^ = ml-m\ = (7.06 - 8.34) x 10"^ eV^ 
lAm^il = \ml - ml\ = (2.13 - 2.88) x 10~^ eVl (2.8) 



The ambiguity in determining the type of neutrino hierarchy is in part due to the large 



amount of mixing in the neutrino sector. As reviewed for example in 25 26 , at the 3a 



level of observation, the magnitude of the entries of the neutrino mixing matrix (2.6) are: 




0.50 -0.63 0.00 - 0.22 
\UpMNs\^{ 0.22 - 0.56 0.44 -0.73 0.57- 0.80 |. (2.9) 

0.40 -0.71 0.59 - 0.82 

Aside from the upper right-hand entry, the content of this mixing matrix has a very different 
structure from the CKM matrix in the quark sector: 



Vckm\ ~ 0.23 0.97 0.04 . (2.10) 




Returning to the parametrization of the mixing matrix given in equation (2.7), the 
current lack of distinguishability between Majorana and Dirac masses implies that there is 
at present no conclusive observational data on the CP violating phases 6,ai and 02- The 
experimental values for the mixing angles have been extracted in |25j, and at the 3a level 
are given by: 

012 ~ 30.5° - 39.3° (2.11) 
~ 34.6° - 53.6° (2.12) 
9i3 ~ 0° - 12.9°. (2.13) 



7 



Current bounds on ^13 from the CHOOZ collaboration |27j are expected to be improved by 
MINOS (28). 

At the la level, global fits to solar and atmospheric oscillation data obtained by KAM- 
LAND and SNO suggest a non-zero value for ^13 29 . In fact, a non-zero value for 613 near 
the current upper bound has recently been announced by MINOS 30 ,^ 



2.3 Neutrinos and UV Physics 

Having described the main experimental constraints, we now review some of the primary 
features of Dirac and Majorana mass terms in the context of the MSSM]^ After this, we 
review the fact that in spite of the suggestive link between neutrinos and high energy 
physics, there is a certain amount of tension in string based models which aim to incorporate 
neutrinos. 

At a theoretical level, there are two features of the neutrino sector which are quite 
distinct from the Standard Model. First, the overall mass scale of the neutrino sector is 
far below the scale of electroweak symmetry breaking, but retains a suggestive link to the 
GUT scale, in that roughly speaking: 

~ (2.14) 

where Auv ~ 10^^ — 10^^ GeV is close to the GUT scale. Second, the mixing angles are far 
larger than their counterparts in the CKM matrix. These observations suggest that neutrino 
Yukawas may have a very different origin from the other couplings of the Standard Model. 

Let us first consider the case of Dirac neutrinos. Simply mimicking the mass terms of 
the Standard Model, the Dirac type interaction: 

W D Xf^HuVN'^ (2.15) 

would then generate a mass term for the neutrinos far above 0.05 eV, unless the entries of 
the corresponding Yukawa matrix are quite small, on the order of 10~^^. This however is 
rather fine-tuned, and it then becomes necessary to explain why all of the other matter fields 
of the MSSM have order one Yukawas, whereas the neutrino sector happens to have such 



small couplings. We will find in section |6j that the relation of equation (2.14) can actually 
be accommodated quite naturally through the presence of a higher dimension operator in 
the MSSM. 



^We thank G. Feldman for bringing this resuh to our attention, which we learned of after the results of 
this paper had already been obtained. 



We refer the interested reader to the review article 31 for further discussion. 
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Leaving aside Dirac neutrinos for the moment, next consider Majorana neutrinos. Al- 



though a Majorana mass term as in (2.1 ) is incompatible with the gauge symmetries of the 



Standard Model, an effective mass term correlated with the vev of if„ can be introduced 
through the higher dimension operator: 



A 



uv 



(2.16) 



where Auv is an energy scale far above the scale of electroweak symmetry breaking. Once 
Hu develops a vev on the order of the weak scale, this will induce a Majorana mass term 



of the type given by (2.1). This operator breaks the accidental global U{1) lepton number 



symmetry of the Standard Model. Assuming that at least one of the eigenvalues of X['^^ is 



an order one number, this will induce the neutrino mass scale of equation (2.14) 



The higher dimension operator of (2.16) can be generated in seesaw models with heavy 



right-handed neutrinos. For example, in the type I seesaw model (considering for simplicity 
the case of a single generation), the superpotential term 



(2.17) 



will induce the requisite effective operator once the heavy Nr field has been integrated out. 
This can be generalized to all three generations of leptons, and to an arbitrary number of 
n right-handed neutrinos (i, j 



1,2,3 and J, J 



1,2, 



,n] 



(2.18) 



While any number of right-handed neutrinos are in principle allowed, in the context of 
four-dimensional SO{10) GUTs the appearance of three copies of Nr is especially natural. 
This is because in addition to the chiral matter of the Standard Model, each spinor 16 of 
SO{10) contains an additional singlet Nr state. Indeed, the presence of three right-handed 
neutrino states renders the U{1)b~l symmetry non-anomalous. However, we note here that 
in the context of string theory, anomalous U{1) symmetries are quite common, and so the 
motivation for precisely three Nr^s is perhaps less obvious. 

While the appearance of a scale close to Mqut is quite suggestive, the bare matter con- 
tent necessary to accommodate the Standard Model and right-handed neutrinos is typically 
insufficient to generate a realistic neutrino sector. For example, although it is a very non- 
trivial and elegant fact that three copies of the spinor 16 in four- dimensional SO{10) GUTs 
contain just the chiral matter of the Standard Model, as well as the right-handed neutri- 
nos, this by itself is not sufficient for generating a Majorana mass term for the right-handed 
neutrinos. Indeed, 16 x 16 is not a gauge invariant operator. 

In four-dimensional 5*0(10) GUT models, it is therefore common to incorporate addi- 
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tional degrees of freedom which can generate an appropriate Majorana mass term for the 
right-handed neutrinos. These extra degrees of freedom can either correspond to additional 
vector-hke pairs in the 16extra © ISextra, or to higher dimensional representations such as 
the 126extra of 5*0(10)]^ The corresponding operators: 

16m X 16m X 16extra X 16cxtra , ^.^ 

W"l6 16l6l6 = 77 , (2.19) 

iWuv 



16 126 — 16m X 16m X 126extra, (2.20) 

can then generate Majorana mass terms for the neutrino component of the spinor once either 
the 16cxtra©16cxtra or the 126extra devclops a vev. In the above, Muv denotes a suppression 
scale which could either correspond to the string or Planck scale. The second possibility is 
quite problematic in the context of string based constructions, since typically, the massless 
mode content will only contain matter in the 10, 16, 16 or 45 of S'O(IO). However, the first 
possibility, involving the presence of higher dimension operators, is compatible with string 
considerations, and has figured prominently in many string based constructions. Note that 
this type of interaction term will also be present in SU{5) GUT models once suitable GUT 
group singlets are included. For a recent example of this type where a suitable combination 
of singlet fields develop vevs, see p2]. 

Even in the context of SU{5) GUT models, selection rules in the effective field theory 
can be quite problematic. For example, in intersecting D-brane configurations, the right- 
handed neutrinos will typically correspond to bifundamentals between two D-brane gauge 
group factors. In such cases, the gauge symmetries of the D-brane configuration forbid 



the coupling NjiNj^. As noted in 33 35 , the additional gauge symmetries of the D-branes 
are often anomalous and so can be violated by stringy instanton effects. Because the 
characteristic size of this instanton is a priori uncorrelated with the size of instantons in 
the GUT brane, an appropriate instanton effect might generate a Majorana mass term in 
the requisite range of 10^^ — 10^^ GeV. Nonetheless, achieving precisely the correct Majorana 
mass scale requires a certain amount of tuning, because the magnitude of the instanton effect 
is quite sensitive to the volume of the cycle which is wrapped by the D-brane instanton. 
Worldsheet instanton effects in compactifications of the heterotic string can also potentially 
generate a suitable Majorana mass term for right-handed neutrinos. 

It is also in principle possible to associate right-handed neutrinos with other GUT group 
singlets, such as moduli fields. In this case, the primary challenge is to obtain a Majorana 
mass which is near the GUT scale. Indeed, moduli stabilization typically will lead either to 
very heavy masses for such fields, or potentially, much lighter masses when one loop factors 
from instanton effects stabilize a given modulus. This is a possibility which does not appear 
to have received much attention in the literature, perhaps because concrete realizations of 

^We recall that the 126 corresponds to the five-index anti-self-dual anti-symmetric tensor of SO{10). 
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the Standard Model with stabihzed moduh are not yet available. 

Even once the correct Majorana mass term has been generated, there is still the further 
issue of addressing more refined features of the neutrino sector, such as mass splittings, 
and the overall structure, or lack thereof, in the neutrino mixing matrix. While it indeed 
appears possible to engineer detailed models of flavor utihzing large discrete symmetries, 
it is not completely clear whether all such features can be incorporated consistently within 
string based constructions. One of the aims of this paper is to show that in a very minimal 
fashion, F-theory GUTs can accommodate mild mass hierarchies and large mixing angles. 

3 Minimal F-theory GUTs 

In this section we briefly review the main features of minimal F-theory GUTs, focusing 
on those aspects of particular relevance for neutrino physics. For further background and 
discussion, see for instance |4(g|9)[lO||T2|[T5)|20], as well as |§]§[Tg[l3||T4||l6}jl9)|2l]. We 
also discuss in greater detail the role of the anomalous global U{1) Peccei-Quinn symmetry 
in the supersymmetry breaking sector of the low energy theory, and its interplay with the 
neutrino sector. 

3.1 Primary Ingredients 

F-theory is defined as a strongly coupled formulation of IIB string theory in which the profile 
of the axio-dilaton tjjb is allowed to vary over the ten-dimensional spacetime. Interpreting 
Tub as the complex structure modulus of an elliptic curve, the vacua of F-theory can 
then be formulated in terms of a twelve-dimensional geometry. Preserving four-dimensional 
= 1 supersymmetry then corresponds to compactifying F-theory on an elliptically fibered 
Calabi-Yau fourfold with a section. In this case, the base of the elliptic fibration corresponds 
to a complex threefold -B3. Within this framework, the primary ingredients correspond to 
seven-branes wrapping complex surfaces in B^. 

In F-theory GUTs, the gauge degrees of freedom of the GUT group propagate in the 
bulk of the seven-brane wrapping a complex surface S, which is defined as a component 
of the discriminant locus of the elliptic fibration. Depending on the type of singular fibers 
over S, the GUT group can correspond to SU (5), or some higher rank GUT group. In this 
paper we shall focus on the minimal case with GUT group SU{5). 

The chiral matter and Higgs fields of the MSSM localize on Riemann surfaces (complex 
curves) in S. The massless modes of the theory are given by the zero modes of these 
six- dimensional fields in the presence of a non-trivial background gauge field configuration 
derived from fluxes on the worldvolumes of the various seven-branes. The Yukawa couplings 
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of the model localize near points of the geometry where at least three such matter curves 
meet 13 

An intriguing feature of F-theory GUTs is that imposing the condition that gravity can 
in principle decouple from the GUT theory imposes severe restrictions on the class of vacua 
suitable for particle physics considerations. This endows the models with a considerable 
amount of predictive power. For example, the existence of a decoupling limit requires that 
the GUT seven-brane must wrap a del Pezzo surface. In particular, the zero mode content 
of the resulting theory does not contain any adjoint- valued chiral superfields, so that for 
example, embeddings of standard four-dimensional GUTs in F-theory cannot be decoupled 
from gravity. Breaking the GUT group requires introducing a non-trivial flux in the U{1)y 
hypercharge direction of the GUT group [5|[8] . The resulting unbroken gauge group in four 
dimensions is then given by SU{3)c x SU{2)i x U{1)y- 

The ubiquitous presence of this flux has important ramifications elsewhere in the model. 
For example, doublet triplet splitting in the Higgs sector can be achieved by requiring that 
this flux pierces the Higgs up and Higgs down curves. In fact, the requirement that the low 
energy should not contain any chiral or even vector-like pairs of exotics also severely limits 
the class of admissible fluxes. 

This rigid structure also extends to the supersymmetry breaking sector. Generating 
an appropriate value for the /i term in F-theory GUTs requires a specific scale of super- 
symmetry breaking \/F ~ 10^ — 10^ GeV, which is incompatible with gravity mediated 
supersymmetry breaking. Instead, F-theory GUTs appear to more naturally accommodate 
minimal gauge mediated supersymmetry breaking scenarios. In fact, the scalar component 
of the same chiral superfield responsible for supersymmetry breaking also develops a vev, 
breaking a global U{1) Peccei-Quinn symmetry at a scale fa ~ 10^^ GeV. The associated 
Goldstone mode then corresponds to the QCD axion. In addition, some of the common 
problems in gravitino cosmology are naturally evaded in F-theory GUTsQ 

As the above discussion should make clear, the framework of F-theory GUTs is surpris- 
ingly rigid. Nevertheless, it is in principle possible to introduce matter content and fields in 
F-theory models to engineer ever more elaborate extensions of the MSSM. Given this range 
of possibilities, we shall focus our attention on vacua with a minimal number of additional 
geometric and field theoretic ingredients required to obtain phenomenologically viable low 
energy physics. 



^As we will explain in subsequent sections, this is only true in the cover theory, before we quotient by 
the geometric action of the Weyl group defined by the geometric singularity. In other words, some of the 
curves may be identified by monodromies, in which case Yukawa couplings can arise at points where only 
two curves meet. See 17 for a recent analysis of such configurations. 

^In 15 , a scenario of leptogenesis in F-theory GUTs based on a non-minimal neutrino sector with 
Majorana masses in the range of 10^^ GeV was studied. We will see later that in minimal implementations 
of F-theory neutrinos, the natural mass scale of neutrinos is somewhat higher. It would be interesting to 
study the associated leptogenesis scenario. 
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It turns out that these minimal ingredients are frequently sufficient for reproducing more 



detailed features of the MSSM. For example, as shown in 12 , minimal realizations of SU (5) 
F-theory GUTs — with the minimal number of curves and interaction points necessary for 
compatibility with the interactions of the MSSM — automatically contain rank one Yukawa 
matrices which receive small corrections due to the presence of the ubiquitous background 
hyperflux. More precisely, the hierarchical structure of the CKM matrix further requires 
the interaction points for the 5// x 10^/ x 10m and 5/^ x 5m x 10m couplings to be 
nearby, suggestive of a higher unification structure. We will revisit this point later when 
we present models with a single i^g point of enhancement which geometrically unifies all of 
the interactions of the MSSM. 

But as noted in [5 , 12 , there are strong reasons to suspect that the neutrino sector of 
F-theory GUTs is qualitatively different. Identifying the right-handed neutrinos in terms 
of modes localized on matter curves, the fact that the right-handed neutrino is a singlet 
of SU (5) implies that the corresponding curve only touches the GUT seven-brane at a few 
distinct points. In |5|, it was shown that Dirac neutrinos could be accommodated from an 
exponential wave function repulsion due to the local curvature of the GUT seven-brane. 
Moreover, it was also shown in [H] that by including additional GUT group singlets which 
develop a suitable vev, it is also possible to accommodate Majorana masses. On the other 
hand, both of these scenarios are somewhat non-minimal in that they require the presence 
of an additional physical input, such as a particular exponential hierarchy in the Dirac case, 
or a new GUT group singlet with a suitable vev in the Majorana case. In this paper we show 
that even without introducing a new scale, or a new set of fields which develop a suitable 
vev, the geometry of F-theory GUTs already naturally contains a phenomenologically viable 
neutrino sector. 



3.1.1 Local Models and Normal Curves 

One of the important advantages of local F-theory GUT models is that some features 
pertaining to Planck scale physics can be deferred to a later stage of analysis. Indeed, 
this is possible precisely because the dynamics of the theory localizes near the subspace 
wrapped by the GUT seven-brane. On the other hand, by including fields such as right- 
handed neutrinos which localize on curves normal to the GUT seven-brane, it may at ffist 
appear that such modes cannot be treated consistently in the context of a local model. As 
we now explain, such normal curves can indeed form part of a well-defined local model. As 
such, they can be consistently decoupled from Planck scale physics. 

Rather than present a general analysis, we discuss an illustrative example. Consider a 
local model of F-theory where the threefold base i?3 is given as an ALE space fibered over 
a base P^. Although the ALE space is non-compact, it contains a number of homologically 
distinct fiber P^'s, which we label as PKw...,P/ n. By, defines a local model with compact 
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surfaces defined by the Pjj-j's fibered over the base P^. The pairwise intersection of two 
such surfaces will occur at a point in the ALE space which is fibered over P^. Identifying 
one such surface as the one wrapped by the GUT seven-brane, it follows that in this local 
model, there are compact curves inside the GUT seven-brane given by a point in the ALE 
space fibered over P^. The model also contains compact normal curves corresponding to 
fiber P(j-)'s which intersect the GUT seven-brane at a point. Hence, modes localized on 
such normal curves can be consistently defined while remaining decoupled from Planck 
scale physics. Although we do not do so here, it would be interesting to study this more 
general class of local models by extending the analysis presented in [2] . 

3.2 U{l)pQ and Neutrinos 

Selection rules in string based constructions can sometimes forbid interaction terms in the 
low energy theory. In the specific context of F-theory GUTs, the U{1)pq symmetry plays an 
especially prominent role in that it forbids a bare /i and Bfi term in the low energy theory. 
Indeed, U{1)pq symmetry breaking and supersymmetry breaking are tightly correlated in 
the deformation away from gauge mediation found in [9]. However, as we now explain, 
the presence of this symmetry can also forbid necessary interaction terms in the neutrino 
sector. After presenting this obstruction, we show that there is in fact a unique alternative 
U{l)pQ compatible with a Majorana scenario. 

3.2.1 Review of Eq and U{1)pq 

An interesting feature of GUTs is the presence of higher rank symmetries. Indeed, these 
symmetries can forbid otherwise problematic interaction terms. For example, in the context 
of the MSSM, it is quite natural to posit the existence of a global U{1)pq symmetry under 
which the Higgs up and Higgs down have respective U{1)pq charges qn^ and qn^- Provided 
that qHu + iHd 0) this forbids the bare /z-term: 



thus providing a partial explanation for why /i can be far smaller than the GUT scale. 
Since the Higgs fields interact with the MSSM superfields, the presence of this symmetry 
then requires that all of the fields of the MSSM are charged under this symmetry. 

In the context of F-theory GUTs, correlating the value of the /i term with supersym- 
metry breaking is achieved through the presence of the higher dimension operator: 



(3.1) 




(3.2) 
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where in the above, X is a chiral superfield which locahzes on a matter curve normal to 
the GUT seven-brane. Here, the X, Hu and Hd curves form a triple intersection and the 
above operator originates from integrating out Kaluza-Klein modes on the curve where X 
localizes. When X develops a supersymmetry breaking vev: 

(X) = X + e'^Fx, (3.3) 

this induces an effective /i term of order: 

/X ~ (3.4) 

As estimated in [o], using the fact that Ayy ^ Mqut, generating a value for the fi term 
near the scale of electroweak symmetry breaking requires \/Fx ~ 10*^ — 10^ GeV [o]. In 
this context, the U{1)pq symmetry can be identified with a linear combination of the 
f/(l) symmetries present on the seven-branes which intersect the GUT seven-brane. This 
necessarily requires that X be charged under U{1)pq with charge: 

Qx = Qh^ + qua- (3.5) 



As explained in [9], this type of structure is quite natural in the context of F-theory 
GUTs and is in fact compatible with an underlying Eq structure. Indeed, decomposing the 
27 and 27 of Eq into irreducible representations of 50(10) x U{1)pq yields: 

D 50(10) X U{l)pQ (3.6) 
27 ^ I4 + 10_2 + I61 (3.7) 
27^ 1_4 + 102 + T6_i. (3.8) 

The MSSM chiral matter transform in the I61, while the Higgs fields transform in the 
10_2- In addition, X transforms in the 1_4. This structure is also compatible with gauge 
mediated supersymmetry breaking, with the messenger fields transforming in the IO2. In 
this context, the U{1)pq charges of the various fields are: 





X 


Y 


Y' 


Hu 


Hd 


10m 




U{1)PQ 


-4 


+2 


+2 


-2 


-2 


+1 


+ 1 



(3.9) 



where in the above, Y and Y' denote the messenger fields of the gauge mediation sector. In 
addition to forbidding a bare /i term, a Z2 subgroup of U{1)pq can naturally be identified 
with matter parity of the MSSM. Indeed, by inspection of the above charges, note that the 
charges of the MSSM chiral matter are all odd, while the Higgs fields are even. 
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The choice of charge assignments obtained by embedding all matter fields in representa- 
tions of Eq is problematic for neutrino models with a Majorana mass term which is induced 
by the F-term (if„L)^/Auv- Indeed, under the charge assignments of line (3.9), this op- 
erator has charge —2. While it is tempting to argue that a suitable vev for the X field 
could generate such a term from a higher dimension operator, note that a non-zero vev for 
X will simply break U{1)pq to the discrete subgroup Z14. Since the operator (if„L)^/Auv 
is not invariant under this discrete subgroup, we conclude that compatibility with a Ma- 
jorana mass term scenario requires incorporating another GUT group singlet with charge 
±2. Once this singlet develops a vev, it is possible to consider models which include this 
higher dimension operator. This is somewhat non-minimal, however, so in keeping with 
the general philosophy espoused in this paper, we shall seek an alternative scenario which 
does not require the presence of an additional symmetry breaking sector, the sole purpose 
of which is to solve a single problem. 



3.2.2 Generalizing U{1)pq 

In a broader context, it is possible to consider more general U{1)pq charge assignments. 
We now show that compatibility with SU{5) GUT structures imposes strong restrictions 
on possible charge assigments. We find that there is essentially a unique alternative U{1) pq 
given by the Abelian factor of SU{b) x ?7(1) C 5*0(10) which is compatible with the re- 
quirements of both supersymmetry breaking and the existence of the operator (iJ„L)^/Auv. 

To establish this result, we begin by asking more generally what U{1)pq charge as- 
signments are compatible with the interaction terms of the MSSM. Assuming that all 
fields in the IQm and 5m have respective U{1)pq charges gio and gg, the interaction terms 
5// X 10m X 10m and 5// x 5m x 10m are compatible with U{1)pq provided: 

1 

qw = -^Qh^ (3.10) 

1 

^5 = -(iHa - Qw = -Qh^ + -^qHu- (3.11) 
If we now demand that the operator (iJ„L)^/Auv is invariant under U{1)pq, we also find: 

+ =0. (3.12) 
Solving for all PQ charge assignments yields: 





X 


Y 


Y' 


Hu 


Hd 


10m 


5Af 


u{iypQ 


+5 


-2 


-3 


+2 


+3 


-1 


-2 



(3.13) 



up to an overall common rescaling of all charges. 
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It is quite remarkable that this structure is in fact compatible with the representation 
theory of SO{10). Indeed, decomposing the 16, 16 and 10 of 5*0(10) into irreducible 
representations of SU{5) x U{1) yields: 

SO{10) D SU{5) X U{1) (3.14) 

16^1_5 + 5+3 + 10_i (3.15) 

16^1+5 + 5-3 + 10+1 (3.16) 

10^5+2 + 5_2. (3.17) 

Note that in contrast to the conventional assignments within the 16 of 50(10), now, 
Hd G 16 and 5m G 10. While this may seem anti-thetical from the perspective of grand 
unification, one of the important features of F-theory GUTs is that locally, the chiral matter 
can organize into the 16, although in the global geometry, this identification is ambiguous. 
Indeed, the localization of interaction terms at points of the geometry can naturally ac- 
commodate both the presence of higher unification structures, as well as the identification 
of this new U{1) symmetry. 

In addition to global symmetries, it is also important to check that matter parity remains 
intact. In fact, with respect to these new charge assignments, note that nothing forbids 
the interaction term HuL. Indeed, because the global U{1)pq symmetry is compatible 
with {HuLy / AuY , it cannot forbid HuL. However, additional discrete symmetries of the 
geometry can in principle forbid such interaction terms, and as proposed in [5], could 
potentially be identified with matter parity. We will return to this point in the context of 
the Kaluza-Klein seesaw in section HI 



3.2.3 F-theory Neutrinos and the LHC 

One of the distinctive features of F-theory GUTs is that integrating out the gauge boson 
of the anomalous U{1)pq gauge theory in general shifts the soft mass terms of the scalars 
away from the value predicted in minimal gauge mediated supersymmetry breaking so that 
the soft scalar mass squared of an MSSM superfield $ obeys the messenger scale relation: 



m$ + AnapQexe^ 



M, 



UWpq 



(3.18) 



where m$ denotes the soft mass in minimal gauge mediation, M[/(i)^q denotes the mass of 
the anomalous U{1)pq gauge boson and apg the associated fine structure constant of the 
gauge theory. In [o] , a particular choice of f/(l) pq charges compatible with an Eq unification 
structure was considered. This choice leads to a predictive deformation away from gauge 
mediation, with potentially measurable consequences at the LHC |9||l0)|20|. Here, we see 
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that considerations from neutrino physics can prefer a different choice of charge assignments 
inducing a different shift in the soft masses. Thus, determining the form of the mass shift 
constrains the form of the neutrino sector, and the converse holds as well! More generally, 
note that we have identified the two Abelian factors in the breaking pattern: 

^6 3 SO{10) X U{l)b D SU{5) X U{l)a X U{l)b (3.19) 

as potential U{1)pq symmetries. However, Majorana neutrino masses single out U{l)a as 
the PQ symmetry. 

Having shown that there is in principle no obstruction to accommodating neutrino 
physics in F-theory GUT scenarios, we now proceed to study the geometry of such config- 
urations. The suggestive link between the suppression scale Auv and the higher dimension 
operator (if„L)^/Auv indicates the presence of GUT scale physics, and so we now turn to 
Majorana neutrinos in F-theory GUTs. 



4 Majorana Neutrinos and the Kaluza-Klein Seesaw 

In this section we study minimal implementations of the Majorana scenario in F-theory 
GUTs. This amounts to determining geometries which contain the terms: 

w.„ 3 ^•i'^qi^. (4.:) 



UV 



As explained in subsection |2.3 this type of operator can naturally originate from a type I 



seesaw mechanism with a superpotential term of the form: 

W D yl'J ■ H^UNi + Mu ■ N^^Ni (4.2) 

for some number of right-handed neutrinos Nj^ labeled by the index /. Here, M/j denotes 
the entries of a Majorana mass matrix. In matrix notation, the coupling A^jV^uv is then 
given by: 

For simplicity, in this section we exclusively consider scenarios where the right-handed neu- 
trinos localize on curves. Indeed, the implementation of bulk mode right-handed neutrinos 
appears less straightforward in the context of F-theory GUTs, although we shall briefly 
comment on this possibility later in the context of Dirac neutrino models. 

Neutrino physics in F-theory GUTs has been discussed previously in |5j. In that context, 
the Yukawa coupling y turned out to be somewhat smaller than an order one number. 
Moreover, upon estimating the expected vev of GUT singlets to be 10^^ GeV, it was argued 
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that an appropriate GUT group singlet P could generate the requisite Majorana mass 
^maj ~ Similar seesaw mechanisms induced by higher dimension operators of flipped 
SU{b) F-theory GUTs have also been studied in [5,36 . 



In a certain sense, however, such scenarios must be viewed as incomplete until we specify 
how P develops a suitable vev. While suitable brane constructions are likely available to 
achieve this goal, in this section we instead investigate minimal constructions which do not 
require additional low energy field theory dynamics. To this end, we demonstrate that in 
F-theory GUTs, it is also quite natural to treat right-handed neutrinos as Kaluza-Klein 
modes. Thus, rather than specify a means by which such fields develop a mass, the fact 
that they are massive modes is already present, by construction. See for example, 37 for 
other scenarios which attempt to realize a seesaw mechanism using heavy modes of the 
compactification. 



In subsection 4.1 we analyze the effective field theory of the Kaluza-Klein seesaw. Since 
the mass term for the Kaluza-Klein modes pairs Nr with N^, while only Nr directly 
couples to the MSSM, we explain why the Kaluza-Klein seesaw requires an identification 
between and N'j^. We then turn to explicit realizations of the Kaluza-Klein seesaw in 
F-theory. We first give in 4.2 a simple realization where the identification comes from the 
geometry itself. We then provide more natural F-theory realizations where the identification 
is provided by monodromy in seven-brane configurations. The main properties of such 



monodromies are reviewed in subsection 4.3 In subsection 4.4 we present a toy model 



based on an SU{7) interaction point which implements the Kaluza-Klein seesaw. This 
example turns out to be only semi-realistic in that it requires the Higgs up and lepton 
doublets to localize on the same matter curve. In subsection 4^ we present a more refined 
example based on an Eg interaction point which accommodates a richer class of interaction 
terms. We also provide a complete model where all interactions of the MSSM unify in Eg. 



4.1 Right-Handed Neutrinos as Kaluza-Klein Modes 

At a conceptual level, it is somewhat ambiguous to interpret the right-handed neutrinos of 
a GUT scale seesaw as zero modes. Indeed, assuming that A^J^ has an order one eigenvalue, 
the resulting Majorana mass scale is quite close to the Kaluza-Klein scale. This observation 
opens up the possibility that right-handed neutrinos may in fact be Kaluza-Klein modes^ 
From this perspective, it becomes unclear whether any right-handed neutrino zero modes 
are in fact necessary. 

Although seemingly quite simple, there is one immediate objection to such a proposal. 
Indeed, when right-handed neutrinos localize on matter curves, the natural expectation is 



^We note that the idea of using Kaluza-Klein modes as right-handed neutrinos has appeared for instance 



in 38 (see also references therein), albeit in a different context. 
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that the six-dimensional field transforms as a bifundamental under the gauge groups of two 
distinct seven-branes. The mass term for the Kaluza-Klein right-handed neutrinos pairs 
Nji and so that the actual interaction term is of the schematic form: 

WDy HuLNr + M^^ ■ NrNI,, (4.4) 

where for simplicity we have included the contribution from a single generation of lepton 
doublets. In this subsection we will suppress all such generational indices to avoid cluttering 



the discussion. To obtain the higher dimension operator of (4.1) after integrating out the 
Kaluza-Klein modes, we would also need a coupling of the form HuLN'^. Note, however, 
that the gauge symmetries of the other seven-branes will forbid such a term! 

The problematic nature of this coupling stems from the fact that the right-handed 
neutrino transforms in a complex representation of the seven-branes. There is, however, 
a more general possibility in F-theory due to the interplay between geometric and field 
theoretic degrees of freedom. In particular, an appropriate discrete group action can identify 
the resulting seven-branes, so that Nr and N'^ transform in a real representation. This 
effect can be analyzed purely in field theoretical terms by passing to a covering theory with 
additional fields which are to be identified in a suitably defined quotient theory. In the 
simplest example, we shall be interested in a covering theory with matter fields L, L', Hu, 
H'^, Nji and N^. The quotient theory is defined by quotienting by the Z2 identification: 

L ^ L', ^ Hi Nr ^ iV^. (4.5) 

Physically, the Z2 group action corresponds to an interchange of the seven-branes under 
which the various bifundamentals are charged. A variant of the seesaw mechanism is present 
in the covering theory provided that Nr couples to Hu and L, with a similar coupling 
between and H'^ and L'„ so that: 

WDy- HuLNr + y'KL'N^ + M^^ ■ iV^iV^ (4.6) 

Integrating out Nr and A''^ then generates an effective superpotential term: 

~ {hJ.){h'J') 

Weff D yy ■ — • (4-7) 

Note that in order for Wefj to be invariant under the exchange symmetry, we must have 
y = y'. Descending to the quotient theory, the corresponding field theory will then contain 
the effective term: 

^eff 3 y VJKK ' (4.8) 
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in the obvious notation. See figure [T] for a depiction of the quiver theory associated with 
these interaction terms in both the covering and quotient theory. This type of mechanism 
clearly extends to multiple generations, and will therefore induce masses for the remaining 
neutrinos. 

It is in fact possible to generalize the field theory construction above to Kaluza-Klein 
seesaw models where we quotient by the action of a more general finite group &, such that 
Nji and A^^^^ are identified. In the covering theory, the action of the finite group & will map 
the matter fields to one another. Let us group the matter fields in the covering theory in 
terms of orbits under the action of &. Consider a covering theory with matter content: 

Orb{H^) = \^a{Hu)W e , (4.9) 

Or6(Z) = |(T(Z)|a G ©} , (4.10) 

OrbiNn) ^ ^^a{N^)\cr e e] , (4.11) 

Orb{N^^) ^ {a{N^^)\a e e] . (4.12) 

In the quotient theory, all fields belonging to the same orbit will be identified. In par- 
ticular, since we want Nji and to be identified by the action of &, we must require 
that Orb{Nfi) = Orb{Nfj). To realize the Kaluza-Klein seesaw, the superpotential of the 
covering theory must contain terms of the form: 

W D y,,, ■ HlL^Nl, + m'fkHiu'iN'kf + M, ■ iV^(iV]j)^ (4.13) 

where the indices i,j,i',j',k label elements in the group orbits such that the associated 
terms form gauge invariant combinations. 

In the context of F-theory, the finite group & will descend from a geometrical symmetry 
of the compactification. Since the matter fields localize on curves, this geometric identifi- 
cation will also identify curves in the covering theory. Thus, distinct matter curves in the 
quotient theory must come from curves in the covering theory lying in distinct orbits of &. 
For example, compatibility with doublet triplet splitting requires the Higgs up and lepton 
doublet to localize on distinct curves. Hence, these fields must lie in distinct orbits in the 
covering theory. 

We now turn to explicit realizations of the Kaluza-Klein seesaw in F-theory. 

4.2 A Geometric Realization of the Kaluza-Klein Seesaw 

In this subsection we present a first realization of the Kaluza-Klein seesaw by directly 
interpreting the ingredients of figure [T] as an intersecting seven-brane configuration which 
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SU {^)gut <SC/(5)gut 

Figure 1: Quiver diagram of the field theory associated to the Kaluza-Klein seesaw in the 
covering theory (left) and the quotient theory (right). 



admits a Z2 group action. Let us first study the covering theory. In terms of the local 
geometry, this can be modelled in terms of two interaction points where SU{5) enhances 
to SU{7) so that the HuLN^ localizes at a point P, while H'^L'N^ localizes at a point P'. 
Geometrically, the required quotienting procedure amounts to the following identification 
of curves: 

so that the Higgs and leptons of the covering theory are correctly identified, while the 
neutrino curve maps to itself. Note that the Z2 group action will in general not leave the 
curve fixed pointwise, even though it is mapped to itself. The intersection points P 
and P' are mapped to each other. 

The covering quiver of figure [T] contains three distinct gauge group factors, which we 
identify with three seven-branes wrapping divisors in the threefold base. Labelling these 
divisors as Tqut-, r+, and r_, the corresponding matter curves are contained in the pairwise 
intersections: 

cr+nrci/T (4.i5) 
S5,,S2, cr.nFct/T (4.16) 
E^cr+nr_. (4.i7) 

Note that while the divisors wrapped by the seven-branes are irreducible, the intersection 
curves may be reducible. Therefore, in this context the lepton and Higgs can in principle 
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SU{5)gvt 




Figure 2: Depiction of a minimal implementation of the Kaluza-Klein seesaw in which the 
fields of the covering theory are identified in the quotient theory. This geometrical action 
also identifies the two interaction points of the covering theory. 



live on distinct curves, and there is a priori no obstruction to realizing doublet triplet 
splitting, although this must be checked explicitly in a given geometric model. The exchange 
symmetry then identifies the divisors and r_, while Tqut remains invariant: 

T+^T_,TGUT^rGUT. (4.18) 

See figure |2] for a depiction of the covering and quotient theory associated with this real- 
ization of the Kaluza-Klein seesaw mechanism. 

As an example which realizes this type of configuration, we consider a configuration of 
colliding A-tjpe singularities. To this end, let u, v, z denote three local coordinates of the 
threefold base such that 2 = is the location of the GUT seven-brane. The z coordinate 
labels the direction normal to the seven-brane, and so the right-handed neutrino curve will 
be parameterized by this coordinate. In this language, the HuLNn interaction point P 
descends from the codimension three enhancement in the singularity type: 

near P : y'^ = x'^ + z^{u — v — a)(u + v — a). (4.19) 

This corresponds to an SU{5) GUT seven-brane at z = 0, and locally defined U{1) factors 
atu — V — a = and u + v — a = 0. Labelling a patch containing P of a given subvariety V 
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as Up{V), the local profile of the divisors, matter curves, and local interaction point are: 



Up (Tout) — 


{z 


= 0), 




Up (r±) = 


{u 


± V — a = 0), 






{z 


= 0)n(M + t'-a = 


0), 


Up (Sj) = 


{z 


= 0)n(M — V — a = 


0), 


{^n) = 


{u 


+ t,-a = 0)n(M- 


V — a = 


Up{P) = 


{z 


= 0) n (w = 0) n (m 


= a). 



0) 



a)n{v = 0), 



(4.20) 



Similarly, the H'^L'N^ interaction point descends from the codimension three enhancement 
in the singularity type: 



near P' : y'^ = x'^ + z^ {u — v — b){u + v — b) , 



(4.21) 



where now we can model the local profile of the divisors, matter curves, and interaction 
point as: 



Up, (Xgut) — 


{z = 0), 




Up, (r±) = 


{u±v — b~0), 




Up, (Sj,) = 


{z^O)n{u + v-b^ 


0), 




(^z ^0)n{u - V - b ^ 


0), 


Up, (s^) = 


{u + v-b = 0)n{u- 


V — 


Up, (P') = 


{z = o)n{v = o)n {u 


= b) 



(4.22) 



Comparing the local data defined by these two patches, we conclude that the Z2 group 
action which interchanges P and P' is given by a reflection in the v coordinate and an 
interchange of the parameters a and b: 



z ^ z, u ^ u 
V —V, a -H- b. 



(4.23) 
(4.24) 



While this provides a simple realization of the Kaluza-Klein seesaw, it is somewhat 
unsatisfactory, in the sense that it requires an additional geometric ingredient to be added 
by hand. Even so, geometries with an appropriate Z2 can in principle be manufactured, 
providing a straightforward realization of the Kaluza-Klein seesaw. Nevertheless, as we 
now explain, there arc other geometric realizations of the Kaluza-Klein seesaw where the 
finite group action comes from monodromy around codimension three singularities. Such 
monodromies occur quite generically in compactifications of F-theory, and therefore provide 
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another means by which to reahze the Kaluza-Klein seesaw. 



4.3 Weyl Groups and Monodromies 

In F-theory compactifications there is a natural set of discrete group actions which are 
especially prevalent, corresponding to monodromies around codimension three singularities. 
In subsequent subsections we will present explicit realizations of the Kaluza-Klein seesaw 
which use the presence of this natural identification. In this context, both the quotient 
and the covering theories will have only one interaction point, in contrast to the geometric 
example of the previous subsection. 

Let us start by reviewing the appearance of monodromies in F-theory seven-brane config- 



urations. We refer the reader to 17 for a recent study of monodromies in compactifications 
of F-theory with codimension three singularities. 

To illustrate the main feature of monodromies, we consider a stack of D7-branes 
wrapping a hypersurface defined by the local equation z = 0. This corresponds to a local 
An_i singularity: 

y^ = x^ + z"". (4.25) 

Assuming that the geometry admits a suitable deformation to lower degree terms, the A^v-i 
can break to Aiy_3 as: 

y^ = x^ + ,^-^(z-h){z-h), (4.26) 
where the ti may be viewed as non-trivial polynomials in the coordinates defined on the 



threefold base B^. Equation (4.26) defines a configuration of {N — 2) D7-branes wrapping 
z = 0, and a single D7-brane wrapping each hypersurface z = ti. In the original SU{N) 
gauge theory, this corresponds to allowing an adjoint- valued chiral superfield develop vevs 
in the Cartan of SU{N) such that SU{N) breaks to SU{N - 2) x f/(l) x f/(l) c SU{N). 
Thus, suitable vevs in the Cartan of the gauge group translate into deformations of the 
corresponding geometry. 



Expanding out equation (4.26), note that we can also write this singularity as: 



= + z'^'-'iz' + az + b), (4.27) 

with a = — ti — t2 and b = tit2- Geometrically, however, it is now immediate that there 
is a broader class of geometries where a and h do not necessarily decompose in terms of 
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polynomial tj's. Indeed, formally solving for the tj in terms of a and b yields: 



*■ = - "^^f^ (4.28) 

= - °-^f^ . (4.29) 

The presence of the branch cut structure indicates that in this more general case, mon- 
odromy around the brane configuration will now interchange the location of the D7-branes 
wrapping z = ti and z = t2. In other words, quotienting by the Z2 symmetry which 
interchanges the two branches, we obtain a single smooth irreducible surface wrapped 
by a D7-brane. This more general breaking pattern corresponds to the decomposition 
SU{N - 2) X SU{2) X f/(l) c SU{N). The branch cut structure reflects a deformation by 
the Cartan subalgebra of SU(2) x f/(l) modulo the Weyl group. This entire system can be 
studied in terms of a covering theory with local coordinates ti and ^2 subject to an overall 
quotient by this Weyl group. Note that we can always study the covering theory, and then 
perform a suitable quotient. 

We now formalize the above procedure of parameterizing deformations in terms of di- 
rections in the Cartan modulo the Weyl group. Starting from a seven-brane with gauge 
group Ggut, we consider an intersection point where the gauge group is enhanced by rank 
r > 2 to G. Let Ggut x G± C G denote a maximal subgroup. The local field theory near 
this intersection point is determined by a theory of deformations of G preserving Ggut- 
Such deformations are parameterized by the Cartan subalgebra of G±, which we denote by 
g±, modulo the Weyl group of G±, which we denote by W{G±) [39] . 

To analyze the action of the Weyl group, we first describe the covering theory where 
deformations are parameterized by g±, and then quotient by W{G±). Let {ti, . . . , t^} 
parameterize directions in g±. Suitable vanishing loci for the tj's then define the locations 
of enhancements in the singularity type of the F-theory compactification. As in [4 ,40 , to 



work out the matter content of the covering theory, we proceed as follows. First, we write 
down the decomposition of the adjoint of G under the maximal subgroup Ggut x G± C G: 

G D Ggut x G± 
adjoint(G') ^ 0(i?i, R[), (4.30) 

i 

where Ri and R[ are irreducible representations of Ggut and G± respectively. Next, to each 
representation R'^ is associated a set of weights, which are points in the dual space g*j_ to 
the Cartan subalgebra. The weights of the representations R[ give the U{lY charges of the 
decomposition of R[ under the branching G± D U^iy to the Cartan subgroup. We have 
thus obtained the matter content in the covering theory. By duality, the weights also give 
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linear combinations of the Cartan parameters {ti, . . . whose vanishing loci define the 
matter curves where matter in the representation Ri of Gqut localizes. The interactions 
arising at such intersection points in the covering theory can be found by writing down 
gauge invariant combinations of the matter content. 

We now study the action of the Weyl group, in order to understand the quotient theory. 
The Weyl group acts on the Cartan subalgebra parameterized by the tj's. In terms of the 
geometry, this corresponds to a group action on the vanishing loci in the geometry, so that 
the Weyl group W{G±) will in general identify some of the matter curves of the covering 
theory. What is particularly interesting is that, as we have just seen, in studying a more 
generic class of deformations of the geometry, such identifications by the Weyl group occur 
generically in compactifications of F-theory! 

In fact, for a generic choice of complex structure all curves of a given Weyl group orbit 
will be identified. In particular, all six-dimcnsional fields with a given representation under 
SU {5) GUT would then be forced to localize on the same matter curve in the quotient theory. 
This turns out to be too constraining for us since, as we noted earlier, the Higgs and leptons 
(which both descend from six-dimensional fields in the 5 © 5 of SU (5) gut) must localize on 
different curves in the quotient theory. However, at the expense of losing a bit of generality, 
it is also possible to consider geometries where only a subgroup of the full Weyl group acts 
to produce the quotient theory. This requires a somewhat more specific choice of complex 
structure in the geometry, since those deformations of the singularity are not fully generic. 

Our discussion of matter curves has been at the level of a quotienting procedure. This 
raises the interesting question to what extent localized modes of the covering theory descend 
to localized modes of the quotient theory. To a certain extent, the notion of localized modes 
depends on the profile of the Kahler metric in both the cover and quotient. As an example, 
consider local coordinates of a covering theory x and y with Kahler form given by: 



Assuming that the quotient acts by sending x —x with y invariant, we now make the 
identifications x^ — x and y = y. The Kahler form of the quotient theory is then of 
canonical form: 



The precise form of the Kahler form in the covering and quotient theories differ, and so 
will affect the extent to which the corresponding modes satisfying the Dirac equation are 
indeed localized along specific loci. Strictly speaking, however, it is not necessary to specify 
the global profile of the Kahler form. Indeed, we shall often be interested in only the 
local profile of modes near a given interaction point. Such effects are controlled by the 
local curvature of the metric and gauge fields, and so we shall typically assume that an 



u! — i(4|x|^ ■ dx A dx + dy A dy). 



(4.31) 




(4.32) 
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appropriate notion of localization is available in such cases. Indeed, in the specific context 
of Majorana neutrino scenarios where we shall consider massive mode excitations anyway, 
the notion of locahzation on a matter curve is itself less well-defined. The important point, 
however, is that an appropriate notion of massive modes with non-vanishing profile near an 
interaction point is still available, and so we shall sometimes abuse terminology and refer 
to "matter curves" in such instances as well. 

Having presented a general discussion of the potential applications of such monodromies 
in seven-brane configurations, we now restrict our attention to some geometric examples. 
As a first toy model, we consider the interaction between L and derived from an 
SU{7) enhancement point. Some deficiencies in this example will then be rectified when 
we present a neutrino sector derived from an enhancement point. 

4.4 SU{7) Toy Model 

As a toy model of the Kaluza-Klein seesaw, we first consider an interaction between the 
Higgs up, lepton doublet and right-handed neutrino curve which originates from a point 
of enhancement to SU (7) in the SU (5) bulk worldvolume theory. Here, the Higgs up and 
lepton doublet localize on two curves where SU{5) enhances to SU{6). The right-handed 
neutrino localizes on a curve which only touches the GUT seven-brane at the SU (7) point 
of enhancement. 

It turns out to be easier to consider instead the parent U{7) theory, locally Higgsed 
down to f/(6) on the matter curves, and f/(5) in the bulk of the GUT seven-brane. The 
analysis is equivalent, but it maintains contact with the perturbative IIB description. 

In the absence of monodromies, the breaking pattern U{7) D U{5) x x U{1)2 

determines three D7-branes with gauge groups U{5), U{l)i and U{1)2 wrapping distinct 
complex surfaces in the threefold base B^. A six-dimensional bifundamental localizes at 
each pairwise intersection of the seven-branes. Two of these bifundamentals localize on 
curves inside of the GUT seven-brane, and may therefore be identified with Hy, and L. The 
final bifundamental is neutral under the U (5) factor and as a GUT group singlet localizes 
on a curve normal to the seven-brane. 

To incorporate the effects of scvcn-brane monodromies, we now pass to a description in 
terms of deformations by the Cartan, modulo the Weyl subgroup. The maximal subgroup of 
U (7) containing U (5) is U{5) xU (2) C U{7). Generic deformations of U (7) preserving U (5) 
are parameterized by the Cartan subalgebra g of U{2), modulo the Weyl group W{U{2)). 
Letting {ei, 62} denote an orthonormal basis, the Cartan subalgebra g is given by the vector 
space {^161 -1-^262}- Thus, the Cartan parameters are {^1,^2}- The Weyl group W{U{2)) is 
the permutation group S2 — Z2 acting on {^1,^2}- 

In this language, the matter content in the covering theory is given as follows. Under 
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the breaking pattern U{7) D t/(5) x U{2), the adjoint of U{7) decomposes as: 

U(7) D U(5) X U(2) 

49 ^(25, 1) + (1, 4) + (5, 2) + (5, 2). (4.33) 

The weights of the vector representations 2 of U(2) are simply and 63 in the dual space 
g*, and the weights of the 2 arc just minus the weights of the 2. Thus, by duality we obtain 
two matter curves in the covering theory where the 5©5 localize, namely ti = and t2 = 0. 
The weights of the adjoint 4 are ±(e^ — 63) and twice the zero weight. There is finally one 
curve where the singlet 1 lives, which is defined by ti = t2- 

Translating into charges under the branching U{7) D t/(5) x ^7(l)i x U{1)2, we obtain 
the decomposition 

U{7) D C/(5) X C/(l)i X U{1)2 

49 ^25o,o + l-i,+i + 1+1,-1 + lo,o + lo,o + 5-i,o + 5o,+i + 5o,-i + 5+i,o, (4.34) 

where the subscripts denote the respective U{1) charges. This recovers the perturbative 
description of three six-dimensional bifundamentals 5_i,o © 5+i,o, 5o,-i © 5o,+i and l+i,-i © 
l-i,+i- 

At the U (7) enhancement point, we find the interaction terms: 

W D 5_i,o X 5o,+i X + 5o,-i X 5+1,0 X l_i,+i + M^^ ■ x (4.35) 

where in addition to the cubic interaction terms derived from the U{7) interaction point, 
we have also included the Kaluza- Klein mass associated with the U{5) singlets. In other 
words, we assume that there are no zero modes transforming in the singlets, but only 
massive Kaluza-Klein modes. Upon making the assignments: 

Hu e 5_i,o, L e 5o,+i, Nr e 1+1,-1, 

K e 5o,-i, L' e 5+1,0, N^R e l-i,+i, (4.36) 
the covering theory superpotential contains the terms: 

W D HuLNr + H'J.'N% + • NrN'^. (4.37) 

Let us now see whether the Kaluza-Klein seesaw mechanism can be implemented in 
this geometric model. To go to the quotient theory, we must quotient by the action of the 
Weyl group. The Z2 Weyl group acts by permuting ti and t2- In terms of U{1) charges, 
it thus permutes the U{l)i and U{1)2 factors. Therefore, quotienting by the Weyl group 
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corresponds to the identification: 



5-1,0 ^ 5o,_i, 5o,+i ^ 5+1^0) 1+1,-1 ^ l-i,+i- 



(4.38) 



Using (4.36), this indeed becomes the required identification 



(4.39) 



It follows that the Kaluza-Klein seesaw will then generate the desired dimension five oper- 
ator in the quotient theory. 

Unfortunately, this toy model is difficult to merge with the other requisite elements of 
F-theory GUTs. The essential problem is that in the covering theory, the fields Hu and 
L' correspond to conjugate representations which localize on the same matter curve, and 
similarly for H!^ and L (see figure [3|. As a consequence, in the quotient theory, and 
L also localize on the same matter curve. This was of course to be expected, since in this 
example we quotiented by the whole Weyl group ^2 = Z2 so that all the curves in the 
covering theory corresponding to the 5 © 5 of SU{5)gut are identified by the Weyl group. 
Doublet triplet splitting of the 5h requires a non-trivial hyperfiux to pierce the Higgs curve. 
This is incompatible with the requirement that the entire 5m SU{5) GUT multiplet is a zero 
mode on the lepton curve. Moreover, if and L' correspond to conjugate representations, 
there is no reason for the bare coupling HuL' to be prevented in the covering theory, since it 
is gauge invariant. This coupling descends to the bare coupling HuL in the quotient theory. 
Therefore, to prevent this bare coupling from appearing in the quotient theory, we must 
again require that Hu and L live on different curves. This geometric sequestering effectively 
plays the role of matter parity. We now turn to an explicit realization of the Kaluza-Klein 
seesaw which incorporates these elements. 

4.5 Eg Kaluza-Klein Seesaw 

In the previous subsection we observed that although generic monodromies in a configura- 
tion of seven-branes would indeed generate a Kaluza-Klein seesaw mechanism, the limited 
number of matter curves in the covering theory always identified the lepton doublet and 
Higgs curve. To overcome this obstruction, we now proceed to configurations with addi- 
tional matter curves in the covering theory. The most straightforward way to implement 
this additional condition is to require that the SU{5) F-theory GUT contain an enhance- 
ment to a higher rank singularity of type G such that the corresponding gauge group 
G D SU{5) X f/(l)^ and Hu, H'^, L and L' localize on four distinct matter curves. 

Adhering to the general requirement that all of the interactions of interest embed con- 
sistently within gauge theory structures, it is therefore most natural to consider the 
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Figure 3: Depiction of the SU{7) toy model described in subsection 4.4 In the covering 
theory (left) Hu and L' localize on the same curve, and the same is true for H'^ and L. As 
a consequence, in the quotient theory (right) Hu and L localize on the same matter curve. 

rank eight singularity G = Eg. Let us now analyze the configuration of curves meeting at 
this Eg intersection point. 

The SU{5) GUT group embeds in E^ in the maximal subgroup SU{5)gut x SU{5)± C 
Eg,. Generic deformations of Eg, preserving SU{5)gut are parameterized by the Cartan 
subalgebra g± of SU{5)±, modulo the Weyl group W{SU{5)±). Let {ei,...e5} be an 
orthonormal basis. The Cartan subalgebra g± of SU{5)± is given by the vector space 
{tiCi + . . . + ^565}, subject to the tracelessness condition X]i=i — 0- These tj's define the 
Cartan parameters. The Weyl group W{SU{5)±) is isomorphic to the symmetric group 
which acts by permutations of the tj's. 

The decomposition of the adjoint representation is then given by: 

^8 dSU{5)gut X SU{5)^ 
248 ^(1, 24) + (24, 1) + (5, 10) + (5, 10) + (10, 5) + (TO, 5). (4.40) 

The Higgs up, lepton doublets and right-handed neutrinos respectively transform in the 5h, 
5m and Iat of SU{5)gut- Hence, they must descend from the irreducible representations 
of SU{5)gut X SU{5)± given as: 

Hue (5,10), Lg (5,10), iV^G (1,24). (4.41) 

Since we will ultimately need to describe the action of the Weyl group of SU {5)± on the 
matter curves of the covering theory, we now identify the corresponding weights of SU{5)± 
associated with each matter curve. Consider the matter curves on which six-dimensional 
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lOcf/T © 10g(7t's of SU{5)gut localize. These fields transform in the 5± © 5^ of SU{5)±. 
Since the weights of the 5_l of SU (5)_l are given by el, 63, Cg, el and el (with opposite signs 
for the weights of the 5_l), it follows that in the covering theory, there are five curves where 
a six-dimensional field in the IOqut ® 10g;7t of SU{5)gut locahzes, specified by the five 
distinct vanishing loci of the tiS. Note, however, that an appropriate choice of fiuxes can 
avoid the presence of any four-dimensional zero modes from such matter curves. 

Next consider matter curves where a six-dimensional field in the 5gut ® ^gut of 
SU{5)gut localize. Such curves will support the Higgs up and lepton fields. These matter 
fields transform in the 10_l © 10j_ of SU{5)±. Since the weights of the 10_l are e* -|- e*, for 
i,j = l,...,5,i^j, there are ten distinct curves with matter content in the 5gut © ^gut 

of SU{5)gut- 

Finally, we consider the right-handed neutrinos. These fields correspond to singlets 
under SU{5)gut, and transform as an adjoint of SU{5)±. The weights of the adjoint are 
e* — e*, i, j = 1, . . . ,5, i ^ j, and four times the zero weight. 

As a result, we obtain that the matter content intersecting at this Eg point lies in the 
following directions of the Cartan subalgebra: 

IOgut ■ ti, i = 1 . . . 5; 

Sgc/t : ti + tj, i,j j; 

Igut- ti-tj, i,j j, (4.42) 

with opposite signs for the conjugate representations. Note that these are subject to the 
tracclessncss condition Yli=iti = 0- For the 1, we omitted the four singlets uncharged 
under U{1)'^. 

This gives us the matter content of the covering theory, where deformations are param- 
eterized by the Cartan subalgebra itself. To each of these directions there is an associated 
matter curve, given by the vanishing locus of the linear combination of the tiS written 
above. Again, we stress that an appropriate choice of fiux can forbid most of these curves 
from acquiring a non-trivial zero mode content. As such, it is enough to focus exclusively 
on the geometric arrangement of curves. In the examples of the next subsections we will 
specify precisely on which curves we allow zero modes. 

The next step is to mod out by the monodromy group to obtain the quotient theory. 
Note, however, that quotienting by the entire Weyl group would identify all the matter 
curves corresponding to a given representation of SU{5)gut- This is not consistent with 
our requirement that the Higgs and lepton localize on distinct curves in the quotient theory. 
Therefore, we will only quotient by a subgroup of the Weyl group; we need to identify which 
subgroup we will be interested in. We first present a simple example where we identify a 
Z2 subgroup of the Weyl group realizing the Kaluza-Klein seesaw with distinct curves in 
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Figure 4: Depiction of the matter curves in the Kaluza-Klein seesaw associated with an Eg 
intersection point. As opposed to the matter curve configuration of figure [3} here and 
L' locahze on different curves in the covering theory (left). In the quotient theory (right), 
Hu and L locahze on two distinct matter curves. 

the quotient theory. After this we present a more involved example in which all of the 
interaction terms of the MSSM unify at the Eg enhancement point. 

4.5.1 A Z2 Model 

Let us first present a Kaluza-Klein seesaw where we quotient by a Z2 subgroup of the Weyl 
group. That is, we consider a geometry where the deformations of the E^ singularity are 
parameterized by the Cartan subalgebra of SU{5)± modulo a Z2 subgroup of the Weyl 
group W{SU{5)±). 

Our aim is now to identify a Z2 subgroup of the Weyl group W{SU (5)_l) which generates 
the Kaluza-Klein seesaw mechanism. We need the following matter fields in the covering 
theory: Hu, H'^ in the 5, and L, L' in the 5, and the singlets Nji and A''^, such that: 

• Hu and H'u lie in a single orbit of the Z2 subgroup; similarly, L and L' form a single 
orbit, as well as Nji and N'^. This ensures that the Z2 subgroup provides the required 



identification of (4.5) 



The two Z2 orbits for the Higgs and lepton doublet are distinct, so that Hu and L 
localize on distinct matter curves in the quotient theory; 

Nji and must have opposite Cartan directions, since they are conjugate fields; 



The Cartan directions allow for the gauge invariant operators given in (4.6) 
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Let us now provide an explicit identification of tliis matter content and Z2 action. We 
consider tlie Z12 subgroup of tfie Weyl group given by tlie permutation (12) (34), wliich acts 
on tlie Cartan parameters as: 



(^1,^2,^35^4,^5) ^ (^2, ^1, ^4, ^3, ^5)- 



(4.43) 



Using the Cartan directions found in (4.42), we make the following matter content assign- 
ments in the covering theory: 



Hu : —ti 
H':-t2 



^4, 



L : t2 + t3, 
L : ti -\- ^4, 



Nr : ti-t2, 
N% : -ti + t2. 



{AAA) 



Note that by construction H^-, H'^, L and L' all localize on different curves, and Nr and 
are conjugate fields. These assignments are consistent with the superpotential terms in 
the covering theory: 



The Z2 permutation (4.43) acts by: 



^K,L^ L', Nr ^ N\ 



(4.45) 



(4.46) 



as required for the Kaluza-Klein seesaw mechanism. See figure |4] for a depiction of this 
interaction structure in the covering and quotient theories, and figure |5] for a depiction of 
how this interaction embeds in a minimal SU (5) F-theory GUT. 

In the above analysis, we have presented one particular choice of Z2 group action and 
matter assigment in the covering theory. In principle, there could be other choices compat- 
ible with the Kaluza-Klein seesaw. In fact, in the next subsection we provide an alternative 
choice which realizes the Kaluza-Klein seesaw and geometrically unifies all the MSSM in- 
teractions at the interaction point. 



4.5.2 Geometric Es Unification of All MSSM Interactions 

In the previous subsection we showed that the Kaluza-Klein seesaw can indeed be acco- 
modated by an Eg enhancement point. In a certain sense, however, it is not particularly 
economical to include such a high rank enhancement simply to incorporate a neutrino sector. 
Indeed, the presence of this higher unification structure suggests that the other interactions 



of the MSSM might also unify at this same point. In fact, as shown in 12 , the hierachical 
structure of the CKM matrix requires the 5h x IOm x 10a/ and Sjj x 5m x 10m interaction 
points to be close to each other. It is therefore quite natural to consider geometries where 
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SU {5)gut 




Figure 5: Depiction of a minimal F-theory GUT with a Majorana neutrino sector. In 
this case, the Higgs up curve forms a triple intersection with the lepton doublet curve and 
the right-handed neutrino curve. Integrating out the massive right-handed neutrino states 
generates the quartic operator (if„L)^/Auv in the low energy effective theory. 
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5'[/(5)gut 




Figure 6: Depiction of a Kaluza-Klein seesaw model in which all of the interaction terms 
geometrically unify at a single point of Eg enhancement. 

all of the interaction terms geometrically unify. 

In this subsection we present a geometry where all MSSM interactions descend from 
a single Eg singularity. In the example we present, only the supersymmetry breaking 
messenger sector localizes at a different point of the geometry. Monodromies play an 
especially prominent role, both in terms of the Kaluza-Klein seesaw, and also through the 
condition that in the quotient theory, the lOjv/'s localize on a single curve. This latter 
condition is important in ensuring that the up type quarks have one heavy generation 
(see [It] for further discussion). To do so, we must however leave the simple realm of Z2 
identifications, and consider the action of a bigger finite subgroup of the Weyl group. See 
figure |6] for a depiction of this geometry with all interaction terms geometrically unified. 

We consider the subgroup & C W{SU{5)±) generated by the order 2 element gi = 
(12) (34) and the order 4 element (72 = (1234). These elements act on the Cartan parameters 
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as: 



91 ■ {tiM,h,U,h) ^ (^2,^1,^4,^3, ^5), 92 ■ (^1,^2,^3, ^4,^5) ^ (^2,^3,^4,^1,^5)- (4.47) 

Note that is invariant under this subgroup. 

We make the following matter assignments in the covering theory. We group the matter 
fields in terms of orbits under the action of ©; we indicate next to each type of field its 
corresponding orbit. 

5h '-{—ti — ts, — 12 — ^4}, 

Sa/ :{^1 + ^2, ^2 + ^3, ^3 + ^4, ^1 + ^4} 

Nr, N% :{±(ti - ta), ±(t2 - h). ±{h - U), ±iU - h)}, 

10m '-{tl, t2,t3, ti}, 

-{tl + ^5, h + ^5, ^3 + ^5, ^4 + ^5}, 

X :{t5 - tl, t5 - t2, h - t3, h - U}, (4.48) 

where = —ti — 12 — — t^^. We include here the singlet X which is required for super- 
symmetry breaking, as explained in section |3| 



Each line in (4.48) corresponds to an orbit under the action of G. Therefore, all of 
the fields in a given line are identified in the quotient theory. By looking at the Cartan 
directions for each field, it is easy to work out all the gauge-invariant interactions in the 
covering theory. The list is rather long, so we will not include it here and instead focus on 
the most salient features. 

First, all Kaluza-Klein neutrino mode and {N'^^ have interactions of the form given 



in line (4.13). Since by construction the 5^, the and the neutrinos live in orbits of the 



finite group, this is sufficient to realize the Kaluza-Klein seesaw mechanism and generate 



the higher-dimension operator (4.8) in the quotient theory. 

Second, after quotienting by the finite group, the MSSM interaction terms x Sa/ x IOa^^ 
and Sjy X IQm x IOa^ are both present in the quotient theory. Moreover, all the IOa^ curves 
in the covering theory are identified by the finite group. Therefore, there is only one IOa^ 
curve in the quotient theory, as required for one heavy up type quark generation. 

Finally, the interaction term XHuHd is not gauge invariant, and so is not present in 
the quotient theory. However, X'^HuHd/Auy is gauge invariant. This operator can be 
produced by integrating out Kaluza-Klein modes, as explained in [9]. Once X develops a 
supersymmetry breaking vev, this generates a /i-term. 

As a result, we obtain an F-theory SU{5) GUT where all the MSSM interactions ge- 
ometrically unify at a single Eg enhancement point. The local effective theory near this 



37 



Eg, interaction point is parameterized by the Cartan of SU{5)±, modulo the particular 



subgroup & of the Weyl group generated by (4.47) 



The only remaining interaction concerns the supersymmetry breaking messenger sector 
jo], which is given by the superpotential term XYY', where Y and Y' are a vector-like pair 
of messenger fields either in the 5 © 5 or in the 10 © 10. Since this interaction has a very 
different origin from the other MSSM interactions, it seems natural not to require that it 
unifies inside the same Eg interaction point. Note however that it could be unified inside 
Eg, although either Y or Y' (or both) would then be required to live on the same matter 
curve as some of the MSSM matter content. 

To end this subsection, we comment that we have not done an exhaustive search for 
finite subgroups of the Weyl group of SU{5)± which allow for the realization of the Kaluza- 
Klein seesaw. There could be other choices compatible with the Kaluza-Klein seesaw. Our 
main goal here was simply to demonstrate that such consistent choices exist. It would be 
interesting to investigate this issue further. 



4.5.3 U{l)pQ and Matter Parity in the Quotient Theory 

Now that we have realized the Kaluza-Klein seesaw, we can analyze symmetries of the low 
energy effective theory directly in the quotient theory. In particular, we can identify the 
U{l)pQ gauge symmetry in the quotient theory, and understand what plays the role of 
matter parity. 

For simplicity, we focus on the two models of the previous subsection with an Eg sin- 
gularity. Recall that in the covering theory, the deformation is specified by the Cartan 
parameters ^1,^2,^3,^4 generating the Cartan subgroup U{iy C SU{5)±. In both models, 
the finite group that we quotiented the covering theory with left the Cartan parameter 
invariant. By the tracelessness condition, we know that ^5 = —ti — t2 — — t^. Therefore, 
all fields in the quotient theory will remain charged under a f/(l) subgroup generated by t^, 
which is the diagonal combination of the four f/(l)'s in the Cartan subgroup of the covering 
theory. It turns out that this invariant U{1) corresponds precisely to the alternative U{1)pq 



presented in subsection (3.2.2) 



Indeed, consider the matter content presented in (4.48) (the same analysis holds for the 
matter content of the Z2 model). The charges of the fields under the diagonal subgroup 
generated by are: 

(4.49) 





X 


Hu 


Hd 


10m 


5m 


Uil)pQ 


+5 


+2 


+3 


-1 


-2 



which are precisely the charges obtained in section (3.2.2). 

We can also say something about matter parity. More precisely, we want to understand 
why the quotient theory admits the quartic superpotential term (if^L)^, while the associ- 
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ated matter parity violating coupling HuL is absent. Note that for the H^L term to be 
present in the quotient theory, one would need terms of the form H^^L^ for some i and j in 
the covering theory. For such a term to be gauge invariant if* would need to have opposite 
Cartan charges to L^; that is, they would need to be conjugate fields living on the same 
matter curve. We may then say that the effective role of matter parity is played by the 
requirement that and L descend from distinct orbits under the action of the finite group. 
It is interesting to note that this requirement is also necessary to implement doublet-triplet 
splitting for the Higgs Hu, which is a priori unrelated to conservation of matter parity. 



5 Yukawas of the Kaluza-Klein Seesaw 

In the previous section we showed that the geometry of F-theory compactifications is flexible 
enough to accommodate a Kaluza-Klein seesaw mechanism, whereby an effective Majorana 
mass for the left-handed neutrinos is induced through a coupling to Kaluza-Klein right- 
handed neutrinos. In this section we estimate the entries of the Yukawa matrix A,-^'' of the 
higher dimension operator: 

Weff^Ki A ' (5-1) 



obtained through the Kaluza-Klein seesaw mechanism. Here, i = 1, 2, 3 is an index for the 
three generations of lepton doublets such that corresponds to the r and doublet. 

This type of interaction term originates from integrating out the heavy right-handed 
neutrinos. In terms of four-dimensional chiral superfields, the Kaluza-Klein seesaw is given 
as: 

W D m,iHj.'Ni + y'^jH^L^N} + MjjN^Nj, (5.2) 

where / and J are indices labelling all of the massive modes of the compactification. In 
matrix notation, the Majorana coupling is then given by: 

^=y.±.^. (5.3) 
Auv M ^ ^ ^ 



The Yukawas yij of line (5.2) are given by overlaps between the Higgs and lepton zero 



mode wave functions with the massive right-handed neutrino zero modes: 

y,j = [ ^HM^^ (5.4) 



Ub 



where Ub C B3 denotes a neighborhood in B3 around the neutrino interaction point. 



The form of the integral in equation (5.4) is to be contrasted with the Yukawas in the 
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quark and charged lepton sectors which are instead given by overlap integrals in a two- 
dimensional neighborhood Us <Z S which contains the corresponding interaction point 
For example, the up type quark Yukawa coupling in the interaction term: 



12 



is given by the overlap integral: 



Us 



(5.5) 



(5.6) 



where the \l/'s denote the corresponding zero mode wave functions. More formally, the 



interaction term of equation (5.6) descends from an appropriate superpotential coupling 



in an eight-dimensional quasi-topological theory. In a perturbative string description, the 



Yukawa of equation (5.4) can be interpreted in terms of holomorphic Chern-Simons theory 



defined in a patch of the neutrino interaction point. We will return to a more precise 



formulation of this overlap integral in subsection (5.2). 



We now explain in crude terms our expectation for the form of the Yukawa matrix y. 
The main point is that whereas zero mode wave functions '^^'^^ satisfy wave equations of 
the schematic form 

= 0, (5.7) 
Kaluza-Klein mode wave functions are massive modes and as such 



(5.^ 



Thus, whereas there is a notion of holomorphicity for zero mode wave functions, there is 



no similar notion for these massive modes. As found in 12 , and as we shall review in 



subsection 5.1, the holomorphicity of the wave function translates into an approximate set 
of f/(l) symmetries which are violated by the presence of background fluxes. These viola- 
tions then generate subleading corrections to the Yukawa matrices of the zero modes. By 
contrast, because the Kaluza-Klein wave functions are not holomorphic, these approximate 
U{1) symmetries will be violated more strongly, leading to milder hierarchies in the neutrino 
sector. 



The rest of this section is organized as follows. In subsection 5.1 we review the computa- 
tion of the Yukawa matrices in the quark and charged lepton sectors obtained in 12 . Next, 



in subsection 5.2, we perform the analogous computation in the case of the Kaluza-Klein 



neutrino sector. With this result in hand, in subsection |5.3[ we compute the form of the 
neutrino Yukawa coupling A(j,) in the low energy effective field theory. Finally, in subsection 



5.4 we discuss the overall mass scale expected from the Kaluza-Klein seesaw, and why the 
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effective seesaw scale can in principle be lower than the GUT scale. 

5.1 Review of Quark and Charged Lepton Yukawas 

As we will explain in the next subsection, the fact that the right-handed neutrinos do not 
correspond to zero modes significantly dilutes the expected mass hierarchy in the neutrino 
sector. To see how this comes about, we first recall the estimate of the quark and charged 
lepton Yukawa matrices obtained in [12] . 

For brevity, we focus on the up type quark Yukawa coupling: 

>if = I '^H.'^hK- (5.9) 

Us 

Although a global description of the wave function profile would be interesting, it is not 
necessary to define the requisite wave functions. Indeed, in a neighborhood of the interaction 
point, the entire gauge theory on S can be modelled in terms of a parent gauge theory which 
is Higgsed down to the bulk gauge group on S by the vev of a locally defined (2, 0) form of the 
parent theory. In the presence of suitable background fluxes, the zero modes correspond to 
solutions to the defect equation of motion of the eight-dimensional quasi-topological theory 
found in [i]: 

uJsAdA^+'-[^,x] =0 + --- , (5.10) 
9aX + ^[0,V'] = O + ---, (5.11) 

where us denotes the Kahler form of the parent theory, denotes the background value 
of the (2, 0) form of the parent eight-dimensional topological theory, x and ip respectively 
denote (2, 0) and (0, 1) forms associated with modes localized on curves where the eigen- 
values of the background (p vanish, and the "■ ■ ■ " correspond to possible higher dimension 
operator contributions induced by background fluxes of the compactification. Here it is 
important to note that the available background fluxes can correspond to local curvatures 
of the metric, gauge field strength, and, for example if-fluxes of the compactification. The 
Yukawa coupling matrix is then defined by evaluating the overlap of solutions to the defect 
equations of motion given above. 

Parameterizing the local patch Us in terms of two coordinates zq and zu such that the 
coordinate of the quark Q curve is zq while that of the U quark curve is zu, the local profile 
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of the zero mode wave functions for Q, U and are [4|[5, 12 : 

- ) ■ exp 













~exp 1 



exp 

exp 

z±z± \ 

' Rl ) 



■ exp (^M 



kl 



ZkZl 



(5.12) 



where we have organized the zero mode wave functions and according to their order 
of vanishing near the mutual interaction point zq = zu = 0. The coordinate corresponds 
to the coordinate normal to the Higgs curve inside of 5, and is given by a linear combination 
of Zq and Zu- Finally, Rq and Ru denote the characteristic lengths of the Q and U curves, 
and R^. ~ denotes the characteristic width of localization, where Mq^j^/M^ 



OtGUT- 



Once we know the form of the wavefunctions, we can estimate the Yukawa coupling (5.9) 



in the presence of background fluxes. The Yukawa coupling is then given by the overlap: 

3-i 



A 



Us 



Rq 



Zu 
Ru 



exp (A^^j ■ ZkZi) ■ Gaussian^ 



(5.13) 



where A^- denotes a quadratic form determined by the background fluxes, and Gaussian 
corresponds to the contribution from the Gaussian wave function factors of the form 
/Rl) such that each six-dimensional field localizes on the appropriate matter 



expi^— \z 
curve. 



In the limit where Aifj is exactly constant, the 3x3 Yukawa matrix is rank one. Indeed, 
note that in this limit the local U{1) x f/(l) rephasing of the coordinates: 



Zq ^ exp{iaQ}ZQ, 
Zu I-)- exp{iau)zu, 



(5.14) 
(5.15) 



causes all Yukawas other than the entry Agg^ to vanish. These f/(l)'s are broken when 
Aifj has non-trivial z dependence. By considering the Taylor expansion of the exponential 



exp (A^jj ■ ZiZj^ , it was proposed in 12 that this "flux distortion" of the wave function 
generates a hierarchical structure in the Yukawa coupling matrix which is characterized by 
the degree of U{1) charge violation^ 



We are now in a position to understand the qualitative difference between the Yukawas 
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^ After 12 and the present paper appeared, much of this hierarchical structure was indeed corroborated 



We shall therefore use the same methodology proposed in 12 to crudely estimate the structure of 



neutrino Yukawas. 
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associated with zero modes and massive modes. As we will see, the internal profile of 
massive modes will always violate this type of U{1) rephasing symmetry. Indeed, massive 
modes are characterized by the fact that (9\E'^^ 7^ 0, so there is no sense in which they will 
preserve the rephasing symmetry present in the zero mode sector. This has the important 
consequence that in computing the overlap between massive modes and zero modes, we 
generically expect to find larger violations of the U{1) rephasing symmetry besides those 
present due to flux distortion. Indeed, since it is a subleading correction, we can safely 
neglect the effects of the flux distortion in the computation that follows. 



5.2 Hierarchy Dilution from Kaluza-Klein Modes 

Having reviewed the estimate of the Yukawa matrices present in the quark and charged 
lepton sector, we now estimate the covering theory Yukawa matrix yij of equation (5.4): 

yi,i = J^hM^^n- (5.16) 

Ub 

To this end, we first discuss the profile of the massive right-handed neutrino excitations, 
and then use this behavior to estimate the form of the neutrino sector Yukawas. 



5.2.1 Massive Mode Wavefunctions 

Since the Kaluza-Klein seesaw crucially relies on the profile of massive modes, we now turn 
to a more explicit description of their internal profiles. There are two ways in which a 
given excitation can correspond to a non-zero mode. The first class of non-zero modes 
descend from massless six- dimensional fields which localize along a given matter curve. 
Upon reducing the profile of these massless six- dimensional field into harmonics of the 
curve, we indeed find massive modes with excitations parallel to the curve. There is, 
however, another class of modes corresponding to fields which are already massive in six- 
dimensions. These turn out to play an especially prominent role in the context of the 
neutrino sector. Since the right-handed neutrino curve is transverse to the GUT seven- 
brane, the corresponding wave functions will have a profile in three directions, corresponding 
to directions transverse to the neutrino curve which we parameterize by the coordinates z± 
and z'j_, and a local coordinate along the neutrino curve. In this section we show that 
the massive six-dimensional fields correspond to harmonic oscillator wave functions in the 
z± and z'j_ coordinates. 

In fact, we can analyze these modes in terms of the quasi-topological eight-dimensional 
theory studied in E] by restricting to the two-complex dimensional patch U given hy z'j_ = 0. 
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The effective action defined over the patch M ' x U contains the terms: 





1^^(0,0) 










/ Tr{ 


^^(0,1) 





2 



t(0'2) (2,0)' 



1 

2 



(5.17) 



where we have included the exphcit Hodge type of each field, and r^'-'^'^^ denotes a zero form 
of the theory. The presence of the neutrino curve is given by the condition that = z^ti, 
where ti denotes an element of the Cartan. Varying with respect to r^^'^''^) and ^/''^'^'^^ we 



obtain the zero mode equations (5.10) and (5.11). 



The background (2, 0) form as well as the bulk gauge field A both play crucial 

roles in defining the zero mode content of the theory. Although the presence of these 
two contributions at first may appear to be on different footings, we note that both can 



be combined in a generalization of the covariant derivative. Indeed, equations (5.10) and 



(5.11) can be written as: 

Va+4>'^ = 0, (5.18) 

where ^ is a vector with entries r]^'^'^\ X^^'^"* ^^id ip^^'^\ and 2^a+<^ denotes the implicitly 
defined differential operator which depends on the background A and (p. In this language, 
the massive modes of the eight-dimensional quasi-topological theory are eigenmodes of the 
Hermitian operator: 

Aa+<p = V\^Pa+^ + Va^^V\^^. (5.19) 

The formal similarity between A and can be made precise using the fact that in 
fiat space, eight-dimensional super Yang-Mills theory originates from the reduction of ten- 
dimensional super Yang-Mills theory. Thus, in a suitably local patch, we can similarly 



view line (5.17) as a reduction of a ten-dimensional theory where simply corresponds to 
a component of the ten-dimensional gauge field. 

Using this observation, we can now deduce the profile of the massive modes in the 
presence of a background A and 0. Since is linear in the coordinate the corresponding 
bulk gauge field in ten-dimensions defines a constant background flux. This leads to a multi- 
dimensional version of Landau's wave function, and as such, the massless and massive 
modes correspond to the ground state, and excited states of a harmonic oscillator in the 
zx direction. Although this corresponds to two real directions, for notational expediency 
we will label the modes in terms of a single integer index > 0. Suppressing all 

dependence on the coordinate zn, the ground state wave function of the harmonic oscillator 
is a Gaussian: 

= exp {-Mlz^Jl) , (5.20) 
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which corresponds to a massless six-dimensional field. The I±-th massive mode is similar, 
and corresponds to exciting the ground state wave function: 

^ f(i^)(z±,zl) exp {-M^z^zl) , (5.21) 

where f^^^^ denotes a degree I± polynomial in z± and 'z±. This corresponds to a massive six- 
dimensional field, with mass set by the characteristic oscillation frequency of the harmonic 
oscillator, so that for I± ^ 0: 

Mj^^o ~ M,. (5.22) 

Note that because contains contributions from two one-dimensional harmonic oscilla- 

tors in the Rez^ and \va.z^ directions, /'^^-^■*(^±, ^) will generically contain contributions 
of all lower degrees as well. 

Returning to the actual case of interest defined by the Kaluza-Klein seesaw, the right- 
handed neutrinos are defined by the vanishing locus z^ = z'^ — ^. Hence, the corresponding 
harmonic oscillator wave functions will now be functions of z^ and z'^. Extending the profile 
of the wave function into the z'^ direction, we thus find that the wave function exhibits the 
profile of a harmonic oscillator in directions transverse to the curve. Letting zi, denote 
the local coordinate along the Icpton curve, which is normal to the neutrino curve, it now 
follows that ^(-'^-l) will contain terms of the form: 

^'-^ D exp {-MlzLZi) (5.23) 

for all i < I±. 



5.2.2 Overlap Between Massive Modes and Zero Modes 

Having estimated the profile of the massive mode wave functions, we now evaluate the 
overlap integral: 

Vi,!^ j^Hjf'M^. (5.24) 

Ub 

Plugging in the rough form of the zero mode profiles for ^ ^"^^ the Yukawa yij is 
then given by: 

Vi,! ~ j d'^ZHd'^ZLd'^ZN (^^^ ' *iv^ i-^ki ■ ^kZi) • Gaussian, (5.25) 

where wc have used the local coordinates for the Higgs, Icpton and neutrino curve zh- 
and zn io define the coordinates of the local patch, and M. denotes the contribution to 
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the profile of tlie wave functions from background fluxes. By inspection, the zh and zn 
integrals give order one answers, up to normalization of the wave functions. Thus, up to 
order one coefficients, the Yukawa yij reduces to an integral over the zl coordinate: 




• exp {Mj^i ■ ZkZi) ■ Gaussian. (5.26) 
Recall that in directions normal to the right-handed neutrino curve, "^^^^ behaves as a 



harmonic oscillator with terms of the form of line (5.23). It now follows that these zl's will 
saturate the overlap integral, and up to normalization of the wave functions, the Yukawa 
yij is then given by: 



\ 3— i 



where Os^i (/) is a step function which is 1 for / > 3 — i and for I < 3 — i. Finally, as 
in [5], the small parameter l/M^,Ri is related to the GUT fine structure constant through 
the relation: 

Writing y as a 3 x matrix where — )■ 00 is the number of massive modes participating 
in the Yukawa, we therefore have: 



(5.29) 




5.3 Neutrino Yukawa Matrix 



In this subsection we estimate the form of the neutrino Yukawas. In matrix notation, this 
amounts to evaluating: 

\H 1 

(5.30) 



A 



uv 



To determine the rough structure of this matrix, consider equation (5.29) in the truncated 
case where y is given by a 3 x A^ matrix with A^ = 4. Summing over all of the massive 
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excitations which have characteristic scale M*, it follows that: 



A^'^) (N = 4) 



1 ^ 



1 



1 




/ 


e 


^1/2 


1 






e 




1 






e 


^1/2 


1 




V 


e 


£1/2 


1 


) 



(5.31) 



(5.32) 



(5.33) 



where each entry of the matrix is multiplied by an order one entry. The generalization to 
an infinite number of modes is now given by 



^(i/) _ 1 ^ L 
Auv ^ M ^ 



e 



£ 

3/2 



,3/2 

e 

a/2 



£ 
,1/2 



(5.34) 



where as before, each entry of the given matrix is multiplied by an order one coefficent. 
Here, the overall coefficient S reflects the normalization due to the contribution of an infinite 
number of modes. In terms of bra-ket notation, the overlap of wave functions leading to 
A^'^V^uv can be written as: 



Auv 



1 



o- 



(5.35) 



In other words, S is specified by the Green's function associated with the massive modes 
of the compactification. 

Finally, we note that the presence of an infinite sum over massive states addresses a 
potential subtlety in that although most of the contributing modes in the infinite sum have 
mass of order M*, a subset of these modes correspond to massless six-dimensional fields. 
These fields descend to massive four-dimensional modes, but with a shghtly lower Kaluza- 
Klein seesaw scale set by the radius of the corresponding matter curve. Because of the 
lower seesaw scale, such modes might at first appear to provide a dominant contribution 
to the seesaw. Note, however, that there is an infinite number of massive six-dimensional 
fields of characteristic mass M*, which overwhelm the contributions from these massless 
six- dimensional modes. 
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5.4 Green's Functions and the Majorana Mass Scale 



In the previous section we obtained a rough estimate for the relative mass ratios in the 
neutrino sector. In this subsection we discuss the overall normalization of the neutrino 
masses set by the heaviest neutrino mass: 

m\ ~ , ^ , (5.36) 

where S is a shorthand for the presence of a suitable regularization scheme defined over 
the infinite modes of the theory. Insofar as M^, is near the GUT scale, the precise value 
of E will determine whether a given geometry will yield a viable mass scale for the light 
neutrinos on the order of 0.05 eV, or will end up being either too large or too small. 

As mentioned previously, S reflects the contribution from the Green's functions associ- 



ated with the massive transverse modes. Returning to the discussion of subsection |5.2.1 
in a patch of the neutrino interaction point, the Kaluza-Klein seesaw can be formulated in 
terms of ten-dimensional fields as: 

W ^ f M'dAM+ f nuCM+ [ KC'U\ (5.37) 

JB-i JB-i JB-i 

where the script fields correspond to ten-dimensional fields, and 9^ denotes the Dolbeault 
operator with respect to the ten-dimensional background gauge field A. Upon reduction 
to eight dimensions, this background corresponds to the background (2, 0) form as well 
as the gauge field of the eight-dimensional quasi-topological theory. Integrating out the 
right-handed neutrinos thus yields: 

f HuC^KC. (5.38) 

Since the neutrino interaction T-LuCM localizes near a point of the threefold base, which we 
denote by P, while 'H'^CM'^ localizes at P', it follows that the resulting term can also be 
written as: 

I HuC5pGa{zb,P')5p'KC' (5.39) 

Jbz 

= H„(P)£(P)G^(P, P')KiP')C'iP'). (5.40) 

More generally, when multiple interaction terms participate in the Kaluza-Klein seesaw, 
the net contribution is of the form: 

L^ff D J]?{„(P)£(P)G^(P, P')K{P')C'{P'). (5.41) 

RP' 
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We conclude that the overall normalization S is 



S~^G^(P,P')- (5.42) 
p,P' 

In particular, when P ^ P', we note that as P — > P', G'^(P, P') diverges, so that: 

P ^ P' : S > 1, (5.43) 

lowering the effective seesaw mass scale. 

We have also seen, however, that in some case the interaction points P and P' coincide 



in the covering theory, as in the enhancement model discussed in subsection 4.5 In this 
case, the geometry near the interaction point P is to be quotiented by the discrete group 
&. Summing over all of the orbits, it follows that in this case, S is given as: 

}^^J2^A{P,cr{P')), (5.44) 

o-ee 

where the limit procedure is defined by taking (y{P') to lie on one of the matter curves 
in the orbit. Note that in this case, the singular behavior of the Green's function will in 
general cancel out, so that in principle, E could be greater or less than one. It would be 
worth investigating this question further. 



6 Dirac Scenario 

Up to this point, we have focussed on Majorana neutrino scenarios. As we now explain, 
the suggestive link between the neutrino, weak and GUT scales is also present in Dirac 
scenarios where the Dirac mass term is generated by the higher dimension operator: 

d'e^S^. (6.1) 

We will show later that this operator is generated in an analogous fashion to the Giudice- 
Masiero operator HuHd/ ^^uy obtained in [9] where Ayy is close to Mqut- Moreover the 
scale of the neutrino mass this leads to is automatically right: Indeed, the most important 
feature of the usual GUT scale seesaw is that: 



m, ~ — , (6.2) 

^^uv ^^uv ^»-uv 
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where as usual, Vu denotes the scale of the Higgs up vev, and Fh^ denotes the F-term 
component of the Hd superfield. Note that Fh^ converts the D-term to a Dirac mass term 
for the neutrinos: 



This last equality follows from the fact that the MSSM superpotential contains the /x-term: 

Wmssm 3 f^HuHd (6.4) 
so that the F-term equation of motion yields: 



dW, 



Fh, - ^^^^^ ~ {Hu) - 10^ GeV^ (6.5) 

where we have used the fact that the /i parameter is typically between 500 — 1000 GeV in 
F-theory GUTs |9]. Similar Planck suppressed operators have been discussed for example 
in 42 . Although the exact operator of line (6.1) was not used, the idea of correlating su- 



persymmetry breaking with the generation of viable Dirac masses has appeared for example 
in |43 



In this section we study minimal F-theory GUT scenarios which incorporate Dirac 
masses through higher dimension operators of the effective theory. As opposed to the case 
of the Kaluza-Klein seesaw, here, the right-handed neutrinos correspond to four-dimensional 
zero modes of the compactification. Moreover, the identification of U{1)pq in S'O(IO) x 
U{l)pQ C -Eg is compatible with Dirac neutrinos. To illustrate the main ideas, we therefore 
restrict to three right-handed neutrino zero modes, and take U{1)pq to be embedded in E^, 
as in [9]. Since the U{1)pq charges of L and Hd are respectively +1 and —2, it follows that 
Np^ has charge —3. Note that the PQ deformation of F-theory GUTs reviewed in section [3] 
will then induce a soft mass term for the right-handed sneutrinos on the order of 100 — 1000 
GeV. 

Although the right-handed neutrinos correspond to zero modes, we will see that Kaluza- 
Klein mode excitations of the higher dimensional theory still play a prominent role in setting 
the overall mass scale of the neutrino sector. 

The rest of this section is organized as follows. In the next subsection we show that 
when the Higgs down, lepton doublet and right-handed neutrino curves meet at a point, the 
required D-term is generated by integrating out massive modes localized on the Higgs down 
curve. In this same subsection we also show that all of the interaction terms of the MSSM 
and the neutrino sector can geometrically unify at a point of E^ enhancement. Next, we 
estimate the form of the Yukawa matrix and find that the resulting mass hierarchy is quite 
similar to the case of the Kaluza-Klein seesaw Majorana scenario. Additional discussion of 
Dirac scenarios in F-theory GUTs is presented in Appendices A, B and C. 
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6.1 Generating Higher Dimensional Operators 



We now demonstrate that the higher dimension operator: 

1 Dirac 



A 



uv 



j d'OHlUN^^ (6.6) 



is generated by integrating out massive modes locahzed on the Higgs down curve. Here, 
the right-handed neutrinos locahze on a curve which is normal to the GUT seven-brane. 
See figure [7] for a depiction of a minimal F-theory GUT which contains a Dirac neutrino 
sector. As we now argue, this operator can originate from an SU{7) point of enhancement 
where the Higgs down, lepton doublet and right-handed neutrino curve form a point of 
triple intersection. 

We begin by writing the relevant terms of the higher-dimensional action in terms of 
an infinite collection of A/" = 1 four-dimensional chiral superfields labeled by points of 
the internal directions of the compactification. As opposed to the Majorana scenario, the 



operator of line (6.6) is obtained by integrating out massive modes localized on the Higgs 
down curve. Treating the higher-dimensional fields as labelled by points of the threefold 
base, the relevant interaction terms are given by: 

ld [ d^e'H\'Hd+ [ d^9n''dCM+ [ d'^en'^dAnd. (6.7) 

The F-term equation of motion for 7{J^ yields: 

dj(Hd + £Ar = 0, (6.8) 

1 

V-d = Hd- —CM, (6.9) 
where Hd denotes the four- dimensional massless mode solution. Plugging Tid back into the 



effective action of line (6.7), we therefore obtain the effective operator: 

\Dirac /■ /" 1 

-f- / d'eH^.UNj, = / d'eHl=rL'Nj,. (6.10) 
^^uv J Jb3 Oa 

In other words, the relevant Yukawa matrix is given by the overlap integral: 

\ Dirac /> -i 

-J—= / ^^^^"^'A (6.11) 
Auv Jbs Oa 

where the \Ef's denote the zero mode wave functions. This can be rewritten in bra-ket 
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SU{5)g\jt 




H\LNr/Kvv 

su{y) 

Figure 7: Depiction of a minimal F-theory GUT with a Dirac neutrino sector. In this case, 
the Higgs down curve forms a triple intersection with the lepton doublet curve and the 
right-handed neutrino curve. Integrating out massive modes localized on the Higgs down 
curve then generates a higher dimension cubic D-term which induces a suitable Dirac mass 
term in the low energy effective theory. 
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notation by inserting a complete basis of states, so that the Dirac Yukawa reduces to a sum 
over massive states I'^n)'- 



Dirac 



(6.12) 

It follows that to estimate the structure of Xf^"^^^ / , it is enough to compute the overlap 
of the massive mode wave functions localized on the Higgs down curve with the lepton 
doublet and neutrino zero mode wave functions: 
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(^Hin^'^) = / ^^n^]., (6.13) 



IUb 

where in the above, Ub denotes a patch in B-^ containing the neutrino interaction point. 



6.1.1 Geometric Unification of All MSSM Interactions 

In the context of the Kaluza-Klein seesaw, we found in section |4] that with an point of 
enhancement it is possible to unify all of the interaction terms of the MSSM at a single 
point of the geometry. In this subsection we show that a similar result also holds for the 
Dirac scenario. See figure |8] for a depiction of this model. 

With notation as in subsection 4.5| we consider a specific discrete subgroup & C 
W {SU{^)±) such that the fields of the covering theory which are to be identified fall into 
orbits of (5. As an explicit example, we consider the order three group generated by the 
cyclic permutation (123) of the first three Cartan parameters: 

(ti, ta, ^3, ^4, ^5) ^ ih, h, ^2, U, h). (6.14) 
In this case, we consider fields in the covering theory which lie in the following orbits: 

5// : {—ti — —^2 — ^3) ~H — ^i} 

^H- {ti+UM + U^h + ti} (6.15) 

W +^5,^2 + ^5,^3 + ^5} (6.16) 

10M:{tiMM, (6.17) 
Nr -.{U - t5}, 

X ■.{U-tiM-t2M-h}- (6.18) 



By inspection, this choice of matter curves allows the interaction terms 5// x 10m x 10 
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We note that in general, (^Hj I'I'-h) 7^ 0. Indeed, this is essentially the content of equation (6.9 1. 
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5'[/(5)gut 




Figure 8: Depiction of a Dirac neutrino scenario in which all of the interaction terms 
geometrically unify at a single point of Eg enhancement. 
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5//x5mx10mi as well as the higher dimension operators X''"if„if(j/Auv and H^LN fij Ayy. 
In this case, the U (1) pq symmetry associated with the abelian factor of 50(10) xU{1)pq C 
Eq is given by the invariant linear combination of [/(l)'s: 

U(l)pQ = C/(l)i + U(l)2 + U(l)s - 3C/(1)4. (6.19) 

One can check that with this identification of t/(l)pQ, we have the charge assignments 





X 


Hu 




lOAf 




Nr 




-1 


-2 


-2 


+ 1 


+ 1 


-3 



(6.20) 



6.2 Neutrino Yukawa Matrix 

In this subsection, we show that the Yukawa matrix of the Dirac scenario has a similar 
structure to that of the Majorana scenario. The zero mode wave functions ^'^ and can 
be organized according to their order of vanishing, so that 

where zl (resp. zn) denotes a local coordinate for the lepton doublet (resp. neutrino) 
curve, and Rl (resp. R^) denotes the characteristic length scale of this curve. As in the 
Majorana scenario, the crucial point is that the massive modes will overlap with the zero 
mode wave functions, inducing maximal violation of the corresponding U{1) coordinate 
rephasing symmetries in the directions transverse to the Higgs down curve. Indeed, the 
massive mode wave function ^^'^^ will contain contributions of the form: 

for all i < II, j ^ In- It now follows that the overlap is given as: 



{^'^''-mK) - V^i"^^"'^3-.(/L)^3-,(/iv), (6.23) 

where here 03-i{I) denotes a step function which is 1 for 7 > 3 — i, and for 7 < 3 — i, and 
we have introduced the small parameters: 

e.^(^]\ e.^(^]\ (6.24) 



Rl J \Rn 
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Summing over all of the massive mode contributions in equation (6.12), it now follows that 
the Dirac matrix is given as: 



\ Dirac 



uv 



( SlEn 
1/2 



1/2 

1/2 1/2 



e 



N 



e 



-1/2 



En 

1/2 
N 

1 



s 



e 



e 

3/2 

e 



3/2 


e 


e 


^1/2 


1/2 


1 



(6.25) 



where S denotes the contribution from the convolution of the wave functions by the Green's 
function, and in the final relation we have used the approximation ^ £■ Comparing 



equations (5.34) and (6.25), we see that the two Yukawa matrices have the same hierarchical 
structure with respect to e. Note that since there is a single neutrino interaction point the 
convolution of the Green's function becomes large near the interaction point. For this 



reason, a similar argument to that given near equation (5.43) of subsection 5.4 implies: 

(6.26) 



S>1, 



which will again boost the value of the Dirac neutrino mass. 



7 Comparison with Experiment 

In previous sections we have seen that with the minimal geometric ingredients required 
to accommodate neutrino physics, both the Majorana and Dirac scenarios yield Yukawa 
couplings in the neutrino sector which are qualitatively different from the case of the quarks 
and charged leptons. In particular, the neutrino Yukawa matrix is given by: 





^3/2 


e 


^3/2 


e 


^1/2 






1 



Ah ~ e^'' e eV2 , (7.1) 
V e 1 / 

where each entry of this matrix is understood to be multiplied by an order one complex 
number. In this section we compare the expected form of these Yukawas with experiment. 
Since we shall mainly be interested in order of magnitude estimates, we will neglect the 
effects of running; this tends to be a sub dominant contribution on top of the theoretical 
uncertainties already present. 

To make contact with experiment, we first extract the expected form of the mixing 
matrix and masses. Depending on the actual geometry of the compactification, the neutrino 
mixing matrix can either exhibit a hierarchy which is milder than that of the CKM matrix, 
or can correspond to a unitary matrix with little hierarchical structure. In both cases, the 
mixing angles are expected to be large so that ^13 should be close to the current experimental 
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bound. Moreover, in the case where the mixing matrix exhibits a hierarchical structure, we 
find the rough relation: 

1 /2 

sin6'i3 ~ sin 6'c ~ chgut ~ 0-2) (7-2) 

where here, 6c denotes the Cabibbo angle. Due to order one ambiguities in the values of the 
underlying parameters, this should be viewed primarily as an order of magnitude estimate. 

After this analysis, we next turn to the expected mass hierarchy in the neutrino sector. 
The neutrino masses exhibit a "normal" hierarchy, with ratios: 

1 /2 

mi : m2 : mg ~ acuT ■ oigut '■ 1- C^-S) 

Again, we find that the milder hierarchy is in reasonable agreement with the observed mass 
splittings. Moreover, using the structure of the Yukawa matrix, we extract the value of the 
lightese neutrino mass mi, and discuss the prospects for testing these expected mass ranges, 
as well as the prospects for distinguishing between the Majorana and Dirac scenarios. 

7.1 Neutrino Mixing Matrix 

As reviewed in section |2| the neutrino mixing matrix is defined by: 

f/PMiV5 = f/f (f/ir^)^ (7.4) 

The Yukawa matrices in the neutrino and charged lepton sectors both exhibit a hierarchical 
structure. However, in terms of the geometry, there is a priori no reason for these hierarchies 
to be manifest in the same basis. Indeed, recall that the hierarchy derives from the presence 
of a local U (1) coordinate rephasing symmetry in a patch of a given interaction point. This 
requires a particular choice of basis for holomorphic functions near this point. Thus, when 
the neutrino and charged lepton interaction points py and pi are far away, there is no reason 
to expect the basis of holomorphic functions to be the same. On the other hand, when py 
and pi are close together, the two basis of holomorphic functions should be approximately 
the same. In particular, if the two interactions occur at the same point, as in the model, 
then the two Yukawa matrices should be in the same basis. This leads to a mixing matrix 
with potentially more structure in the parameter e. 

Thus, the neutrino mixing matrix depends on whether the two interaction points are 
nearby or far away in the geometry. In the following subsections we further discuss these 
two possibilities. 
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7.1.1 Hierarchical Mixing 



As we alluded to previously, the form of the CKM matrix found in 12 strongly hints at 
the presence of a higher unification structure. As noted in [1^, a hierarchical structure in 
the CKM matrix requires the up and down type interaction points pu and pd to roughly 
satisfy the relation \pu — Pd\ ^ 0.1 x M^^j,. Turning the discussion around, the hierarchy in 
this sector can then be taken as evidence for the existence of a higher unification structure. 
Unifying neutrinos with the remaining matter content of the MSSM, it is then natural to 
perform the further identification Pu = Pd = Pi = Pu- Indeed, for both the Majorana and 
Dirac scenarios we presented models of this type, where all of the interaction terms unified 
in a single point of enhancement to E^. 

When pi is close or equal to p^, the rephasing symmetry of the local coordinates will 
lead to a hierarchical structure in the neutrino and charged lepton Yukawas with respect 
to the same basis. These Yukawas are then estimated to be: 




where the form of A(i) was found in 12 . Introducing matrices Ul and f//j such that 
^i'^''A(i/) (ul^^^ and uf^\{i) {u^^^ are diagonal, we note that since a matrix with entries 
\ij ~ has (f/i).. ~ {Ur)^. ~ el^'-'^il, we obtain: 




1 £2 ^4 
r(0 „, I .2 1 ^2 



Ul> ~ £2 1 ^2 1 (7_Q) 




The resulting form of the PMNS matrix is then dominated by the terms in u'f'^ so that: 



U'p^rNs = U^^{ui^)' -\ 1 \ (7.7) 



Here, we have simply estimated each matrix element to be an order one complex number 
multiplied by the appropriate power of e. The diagonal entries of Up^^g are expected 
to be order one complex numbers, so that in the limit where e — )■ 0, Up^^g tends to a 
diagonal unitary matrix. To be more precise, since the off-diagonal entries are small but not 
infinitesimally so, a more reliable estimate of the diagonal entries is obtained by imposing 
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the constraint that Upj^/j^g is unitary, which imphes that: 



1=1 



(7. 



for I = e, fi, T. 

In the context of F-theory GUTs, the parameter e 
this value, we obtain the final expected form for the neutrino mixing matrix: 



M^ut/M^ ~ aj^^. Plugging in 



jjF-th 
'-^PMNS ~ 



a 



\ 



a 



1/4 
GUT 
1/2 
GUT 



1/4 
^GUT 

1/4 



a 



GUT 



O^GUT ^ 



a 



1/2 
GUI 
1/4 



GUT 
Ur3 J 



(7.9) 



where the f/'s along the diagonal are fixed by (7.8). More precisely, (7.9) provides an 
estimate for the magnitudes of the entries of the neutrino mixing matrix. 

It is interesting to compare this form of the neutrino mixing matrix with that of the 
CKM matrix obtained in (121: 



•t/F-th 
^CKM ~ 
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1/2 
^GUT 
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1/2 
^GUT 
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£2 






K 


3/2 

OLQljrp 


acuT 



3/2 \ 



^GUT 
OlGUT 
1 



(7.10) 



which is manifestly more hierarchical. 



1/2 

Let us now compare with experiments. We know that a^W ~ Plugging this value 



into our estimate for the neutrino mixing matrix (7.9) and extracting the diagonal t/'s 



using the unitarity constraint (7.8), we obtain the rough estimate for the magnitudes of the 
mixing matrix elements: 



IrrF-th I 
\'-^PMNs\ ~ 



^0.87 0.45 
0.45 0.77 
0.2 0.45 




(7.11) 



This is to be compared with the experimental result which was quoted in section [2] 



I U^'^ I 



0.77-0.86 
0.22 - 0.56 
0.21-0.55 



0.50 
0.44 
0.40 



0.63 
0.73 
0.71 



0.00 
0.57 
0.59 




(7.12) 



These two matrices look amazingly close! Given that we are working only up to order 
one coefficients, this reveals a very interesting match between theory and experiment. In 
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fact, as in the mixing matrix seems to be relatively insensitive to the various order one 
coefficients which appear in the Yukawa matrices, since the results appear very close to the 
actual experimental result. It would be interesting to study in more precise terms whether 
these order one effects indeed tend to cancel out. 



Given the rough numerical values of the mixing matrix in equation (7.11), we can also 
extract estimates for the values of the neutrino mixing angles. Here it is important to stress 
that the theoretical uncertainties present will mean that the numerical values of the angles 
thus obtained should only be treated as crude approximations. To start with, from the 
form of the matrix we see that 612 and 623 should take similar values, while ^13 should be 
smaller. We roughly estimate 

^13 ~ aJuT -^c- 0.2, (7.13) 



where 9c is the Cabibbo angle, the value of which we have extracted from (7.10 ). Converting 
from radians to degrees, this yields the rough expectation 6'^^"*'^ ~ 10°, where we have 
rounded to the ffist significant figure since order one uncertainties in the coefficients of the 
mixing matrix elements will also propagate to the mixing angles. Since the experimental 
upper bound on ^13 is on the order of 13°, we conclude that from F-theory we expect 6*13 
to be close to its upper boundp] 

We can also extract values for the two other mixing angles. From the entries 12 and 23 



of the mixing matrix (7.11) we obtain *^ ~ ^23 *^ ~ 30°, where we have again rounded 



to the nearest significant figure. These order of magnitude estimates are to be compared 



with the experimental values extracted in 25 which at the 3ct level are: 9i2 ~ 30.5° — 39.3° 
and 923 ~ 34.6° - 53.6°. 

7.1.2 Non-Hierarchical Mixing 

Although somewhat counter to the notion of unification, it is in principle also possible to 
consider geometries where the neutrino and charged lepton interaction points are not close 
together. In this class of geometries, the matrices f/^ and U^^ are hierarchical, but in 
different bases. As a consequence, our expectation is that the mixing matrix Upmns should 
consist of a "generic" unitary matrix with no particular structure. Although we do not have 
a precise notion of genericity, as a substitute we can consider Upmns to be a random unitary 
matrix. To generate random unitary N x N matrices, one has to use the only probability 
measure which is invariant under U{N) group multiplication, known as the Haar measure. 
Perhaps surprisingly, even in this case where little structure is present, we still obtain the 



^"'^ After the results of this paper had already been obtained, we learned from G. Feldman that current 
results from the MINOS collaboration have indeed found evidence that the mixing angle 6*13 is non-zero, 



and is close to this upper bound 30 
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Figure 9: Probability density functions for the three mixing angles 612, 023 and ^13, for a 
random unitary neutrino mixing matrix. 

qualitative expectation that the mixing angles 612 and 623 should be comparable, while 6*13 
should be somewhat smaller. 

This directly follows from the parametrization of the neutrino matrix in terms of the 
mixing angles 6ij. It is at first tempting to think that generating uniformly distributed 
random mixing angles and CP violating phases will generate random unitary matrices 



through the parametrization (2.7). This is however too naive. In Appendix D we review 
the parametrization of the Haar measure in terms of the three neutrino mixing angles. With 
respect to this measure, the probability density functions for the mixing angles are given 



by equation (D.6) of Appendix D: 



P(^^i2) = 2sin(^i2)cos(^i2), 

P(e23) =2sin(^23)cOs(^23), 

P(^^i3) =4sin(^^i3)cos3(^^i3). (7.14) 

This means that to generate random unitary matrices in terms of mixing angles, we should 
not consider uniformly distributed mixing angles, but rather the probability densities of 



equation (7.14). This may seem surprising at first, but is again simply a consequence of 
the way that the neutrino mixing angles parameterize a unitary matrix. See figure [9] for a 
plot of the probability densities for the three neutrino mixing angles. 

From the explicit form of these densities, we extract that the two mixing angles 612 
and 6'23 behave similarly, while the mixing angle 613 has a drastically different probability 
density. In fact, from figure |9] one can see that both distributions of 6*12 and 6*23 have mean 
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value 45°, while the distribution of ^13 has a lower mean value at 33.75°. Therefore, given a 
random neutrino mixing matrix, we expect that 6^ and 623 should be roughly of the same 
order, while 613 should be smaller. This fits relatively well with the current experimental 
data for the mixing angles reviewed in section |2j Amazingly, randomness itself provides an 
explanation why 613 should be smaller than the two other mixing angles! 

However, the actual experimental values are somewhat smaller than the mean values of 
the probability distributions. It is therefore worth asking what is the probability that the 
angles have their measured values, using the probability distributions relevant for random 
unitary matrices. From simple integration of the probability densities shown in figure |9| we 
obtain the following probabilities: 

P(^i3 < 13°) = 9.9%. (7.15) 

We conclude that randomness of the neutrino mixing matrix is potentially consistent with 
the experimental values, provided that 6^13 is close to its current upper bound. For instance, 
if the upper bound was lowered to 1°, we would get the probability: 

P(^,3 < 1°) = 0.06%, (7.16) 

which illustrates the general point that we expect ^13 to be as close to the current exper- 



imental bound as possible. In figure [10} we provide a plot of the probability that ^13 be 
lower than a given angle (the cumulative distribution function), and compare with the same 
probability for the other mixing angles. 

We can also say something about the CP violating phases. As explained in Appendix D, 
the phases of a random unitary matrix are uniformly distributed between and 27r. These 
correspond to the CP violating phases 6, ai and a2 of the neutrino mixing matrix, where 
the latter two are physical only in the Majorana case. Hence, for this class of geometries in 
F-theory there is no reason why these phases should vanish; all values are equally probable. 

As a cautionary note this analysis should only be viewed as semi-quantitative. This 
is because in a completely specified compactification, the neutrino mixing matrix is not 
exactly random, since the actual overlap of all wave functions can be computed. Thus, in 
a more precise computation randomness should be supplanted by an appropriate notion of 
genericity. Even so, this analysis provides a first estimate of what to expect. 
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Figure 10: Plot of the probability that a random unitary 3x3 matrix has a value of 9ij 
less than a given cutoff 9, where < 9 < 90° (plotted up to 30°). The dashed vertical line 
denotes the current experimental bound of roughly 13°. 

7.2 Neutrino Mass Hierarchy 

We now turn to the mass hierarchy in F-theory GUTs. Returning to the form of the Yukawa 
matrix, the neutrino masses are: 

Mgut 

where S is the contribution from the Green's function near the interaction point. Here, we 
have absorbed the overall volume normalization from the Higgs up and lepton doublet wave 
functions into the value of S. These normalization factors tend to decrease the couplings 
by an additional factor of 0{aGUT)- Let us now estimate the overall normalization of the 
heaviest neutrino. Since the neutrinos exhibit a normal mass hierarchy, we have that 



Am|^ = y m\ — mf ~ ~ 50 meV. (7-18) 
For the Majorana and Dirac scenarios we have the two slightly different estimates 

Mgut 
Mgut 
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Using the values Vu ~ 170 GeV, Mgut ~ 3 x 10^^ GeV, and ~ 500 - 1000 GeV |9|, 
it follows that the overall enhancement from the Green's function and volume dependence 
should be SMaj ~ 50 and SDirac ~ 10. This seems to slightly favor the Dirac scenario. 

Having discussed the overall mass scales, we now turn to the relative neutrino mass 
ratios, which from lightest to heaviest are: 



mi : 1712 '■ TTT'S ~ £^ : e : 1. (7-21) 
Compare this with the parametric e dependence in the masses of the charged leptons and 



quarks obtained in 12 : 



rrif. : : ~ : £^ : 1, 
ruu : rric ■ rut ^ : : 1, 
rud : mg : nii, ~ : : 1. 



(7.22) 
(7.23) 
(7.24) 



It is clear that the neutrino hierarchy is much milder than that of the quark and lepton 
sectors. In fact, the neutrino mass hierarchy turns out to be the fourth root of the hierarchy 
in the charged lepton sector sector. In terms of numerical values, recall that the parameter 
£ is roughly estimated to be 

0.2. (7.25) 



1/2 



We note however that a more precise estimate of e will depend on details of the geometry. In 
particular, as explained in [l2] , the overlaps of distinct wave functions will lead to different 



possible values of e. In other words, in equations (7.21 )-(7.24), the parameter e may in 
principle take different values. 

The mass splittings are given by: 



Am: 



2 2 2/1 J 

31 — mg — ~ ""^3(1 ~ "31 ■ ^ 



Am2i 



21 ■ £ 



(7.26) 
(7.27) 



where the (i's are order one coefficients. The mass splittings then satisfy the relation: 



Amli 
Amii 



(7.28) 



Using the values for the mass splittings reviewed for example in 25 , 26 , the mass splittings 
at the 3a level are: 





Min 


Central 


Max 


AmlJilO-^ eV^) 


7.06 


7.67 


8.34 


Am2,/(10-3 eV2) 


2.13 


2.49 


2.88 



(7.29) 
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The maximum and minimum values of equation (7.28) compatible with this range of ex- 
perimental values are then given as: 



min Am2i 
max Amg]^ 



< 



< 



max Amn 



21 



min Am|]^ 



(7.30) 



Extracting the minimal, maximal and central value of e consistent with equation (7.28) we 
therefore obtain the following range of experimental values for e: 





Min 


Central 


Max 


e 


0.16 


0.18 


0.20 



(7.31) 



We note that this value derived from experimental results is consistent with the rough 

1 /2 

estimate e ~ (^gut ~ O-^' Taking into account that our estimates are only valid up to order 
one numbers, this reveals a very interesting match with experiments. 

We can now use the value just obtained to estimate the lightest neutrino mass mi. For 
simplicity, we use the central value e ~ 0.18 obtained in (7.31). Given our prediction that 
neutrino masses exhibit a normal hierarchy we can take and m2 given by the square 
root of the two mass splittings: 



^31 



50 ± 4 meV 
8.7 ±0.4 meV. 



(7.32) 
(7.33) 



F-theory. We obtain 



12 



Using the relative mass ratios (|7.21|), we can then extract the value for nii predicted by 

(7.34) 



3 meV. 



Hence, we expect a relatively light value for mi. As before, we note that this should be 
interpreted only as a crude estimate on the value of mi, since we are working up to order 
one numbers. The expected deviation in the value of e was estimated by comparing with 
the variation present, for example, in the charged lepton sector. Fitting the masses of 



the electron, muon and tau to the hierarchy present of equation (7.22), the central value 
of El ~ 0.36. By contrast, using just the mass ratios m^/rrir ~ ej^ and rrie/m^ ~ e\ 
respectively yield e ~ 0.49 and e ~ 0.26. Thus, we can expect an error of roughly 50% in 
extracting the value of e in the neutrino sector. 



^^It is interesting to compare this value to the landscape of AdS^ Standard Model vacua, where it was 
found that in a Dirac scenario, an AdS^ minimum requires a lightest neutrino of mass greater than 8.3 
meV, a metastable dS^ minimum requires a mass of 7.1 — 8.3 meV, and no minimum is present for lighter 



masses 44 
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7.3 Distinguishing Majorana and Dirac: Neutrinoless Double 
Beta Decay 

In this paper we have found that both the Majorana and Dirac neutrino scenarios naturally 
fit within the F-theory GUT framework. Indeed, up to multiplicative phases which cannot 
be removed in the Majorana case, the form of the Yukawas are identical. It is in principle 
possible to distinguish between these two scenarios through neutrinoless double beta decay 
experiments. Here, the idea is that a nucleus with Z nucleons and atomic number A can 
undergo two decays. The emitted neutrinos can then annihilate each other provided a 
Majorana mass term couples these two states. The associated decay rate for neutrinoless 
double /3 decay scales with the Majorana mass as: 

r ((Z, A) ^{Z,A + 2)e-e-) ~ M • ■ \mpp\^ , (7.35) 

where Anuc is the contribution to the amplitude from nuclear processes, M is an overall 
phase space factor and the relevant Majorana mass term is given by: 



(7.36) 



As reviewed in 45 , experiments such as CUORE 46 , GERDA 47 and Majorana 48 , 49 



will hkely be sensitive to this decay rate at the level of mpp ~ 50 meV. To give a rough 



expectation for future experiments, as explained for example in 50 , within ten years the 
EXO experiment is expected to be sensitive down to nipp ~ 4 — 40 meV. 

In neutrino models with a normal hierarchy, it is sometimes common to also posit that 
the mixing angle 6*13 is quite small. In such models, the value of would instead 

be controlled by nii and m2, rendering this effect less observable. For this reason, it is 
common to say that observing neutrinoless double beta decay would appear to favor a 
Majorana scenario with an inverted hierarchy. 

But in the context of F-theory GUTs, we have seen that it is natural to expect as large 
a value of Oi^ as possible, consistent with current experimental bounds. It is therefore 
of interest to study whether upcoming neutrinoless double beta decay experiments will be 
sensitive to the value of expected in the class of models studied in this paper. Including 
the effects of the CP violating phases, we obtain: 



Depending on the relative phases of these contributions, the individual summands can 
either add constructively, or destructively. For concreteness, we take the rough numerical 
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estimate for the magnitudes of the mixing matrix elements obtained in equation (7.11 ), with 



the values of the masses given in equations (7.33) and (7.34), and range over the values of 



the CP violating phases. As a function of ai, a2 and 5, our rough estimate for mp,^ is then: 

\ml-^'^\ = I (1.4 ± 0.8) -e^"! + (2.8 ± 0.8) ■ e*"^ + (1.3 ± 0.9) ■ e'^^^l meV. (7.38) 

Thus, the maximal value of mpfs expected is: 

m;^;" ~ 6 meV, (7.39) 

while the minimal value expected is consistent with zero. For generic complex phases, we 
therefore roughly expect m^p on the order of a few meV which is too small for observation in 
the current round of experiments, but which is tantalizingly close to the limits of sensitivity 
expected in the near future. 

7.4 Single Beta Decay 

Although the lightest neutrino mass mi we have found is likely to be too small for direct 
detection, it is still of interest to consider constraints from other experiments. Here we 
focus on constraints derived by precisely measuring the masses of all of the visible decay 
products in single beta decays. The effective mass of the electron neutrino, or m/j is: 



\m,f = J2m^U[,^^'\\ (7.40) 



Results from the Troitsk experiment 51 and Mainz neutrino mass experiment 52 give only 



an upper bound of 2.5 eV and 2.3 eV, respectively. The KATRIN experiment is expected 



to be sensitive to a non-zero value of m^j down to 0.2 eV 53 



Again using the rough numerical estimate for the magnitudes of the mixing matrix 



elements obtained in equation (7.11 ), with the values of the masses given in equation (7.34), 
we obtain: 

|mj-*''| ~ 5 - 10 meV. (7.41) 
which is far too small to be observed by current direct detection experiments. 



8 Conclusions 

In this paper we have studied minimal realizations of F-theory GUTs which contain a 
neutrino sector. We have found that small Majorana and Dirac neutrino masses can be 
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accomodated naturally in minimal F-theory GUTs. In both scenarios, Kaluza-Klein modes 
play a prominent role. Owing to the fact that these massive Kaluza-Klein wave functions 
are not holomorphic, the neutrino Yukawa matrix exhibits a milder hierarchical structure 
than its quark and charged lepton counterparts. 

For both the Majorana and Dirac scenarios, we have found that a normal hierarchy is 

1/2 

expected, with relative mass ratios mi : m2 : ~ ogut '■ (^gut ■ which is consistent 
with experimental values for the neutrino mass splittings. In analyzing neutrino mixing, we 
have considered geometries where the neutrino and charged leptons unify at a single point, 
as well as configurations where these interactions do not unify. In the former case, we find 
a mild hierarchy in the mixing matrix compatible with current observational constraints. 

In particular, we find that the mixing angle 613 and Cabibbo angle are related through 

1/2 

^13 ^ 6c ^ '^GUT ~ 0.2. When the neutrino and charged lepton interactions do not unify, 
we instead expect a generic neutrino mixing matrix. For this reason, such models naturally 
realize large mixing angles. This in turn leads to the expectation that the mixing angle 
6'i3 is close to the current experimental bound. In the remainder of this section we discuss 
further potential avenues of investigation. 

In our implementation of the Majorana scenario in F-theory GUTs, it was necessary to 
consider an alternative choice of Peccei-Quinn symmetry consistent with the Majorana mass 
term. It would be interesting to investigate in more detail the phenomenology associated 
with this new choice of Peccei-Quinn symmetry, and in particular the expected form of 
LHC signals, much as in the analysis of [20 j. 

In this paper we have also presented Majorana and Dirac scenarios where all of the in- 
teractions of the MSSM geometrically unify at a single Eg point of enhancement. Our main 
purpose was to present examples in which the monodromy group appropriately identifies 
curves in the quotient theory. Studying other subgroups of the Weyl group of SU (5)_l in the 
breaking pattern SU{5)gut x SU{5)± C Eg may provide further insight into reahzations 
of F-theory neutrinos. 

Finally, in extracting the neutrino Yukawa matrices, the overlap between massive modes 
and zero modes enters in a crucial way in both the Majorana and Dirac scenarios. It would 
be worth studying more precisely how these massive modes fit into the eight-dimensional 
quasi-topological field theory framework, and how the general equations of motion can be 
deformed to include massive excitations. 
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Appendices 

A Dirac Scenario Operator Analysis 

In this Appendix we discuss in more general terms cubic and quartic operators which can 
potentially generate a viable Dirac neutrino mass term. As in section [6} we restrict attention 
to U{l)pQ charge assignments compatible with the embedding 5*0(10) x U{1)pq C Eq. 

At a minimal level, generating a Dirac mass for neutrinos requires the presence of an 
operator which contains the product L^N^^ for i, j = 1, 2, 3, where the N^^ are right-handed 
neutrinos, as well as some additional fields which develop a suitable vev to generate a Dirac 
mass term. Compatibility with SU{2)l x U{1)y gauge invariance therefore requires L^Nj^ 
to couple to either Hu, or iJ^. In keeping with the requirements of a minimal matter 
spectrum, the only fields which develop a vev are H^, and X, where this last field 
develops a supersymmetry breaking vev 

{X)=x + e^Fx (A.l) 

with X ~ 10^^ GeV and Fx ~ 10^^ GeV^. For this reason, we shall restrict our attention 
to operators containing the fields L, Nr-, Hu-, Hd and X. 

At the level of cubic terms in superfields, the possible invariant terms are: 

Oh.iw, = j A^OH.VN'^, (A.2) 



Using the U{1)pq charge assignments described in subsection 3.2.1 obtained from iden- 
tifying U{l)pQ as the abelian factor of S'O(IO) x f/(l) c Eq, it follows that must have 
charges: 

j d^OHuLNR =^ qpQ (Nr) = +1, (A.4) 
d^^^^f^ ^ (^^) ^ _3_ (A.5) 

In other words, in the first case Nr comes from the 27 of Eq whereas in the second case 
Nr comes from the 78. 
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Assuming that the overall coefficient of each operator is an order one number, note that 
the resulting Dirac mass in each case is: 



170 GeV , (A.6) 



rriH^LN, oc ^ ~ ~ 0.01 eV, (A.7) 

where in the second line, we have set Auv ~ 10^® GeV, and used the value of the /i parameter 
obtained in jo] so that: 

OWmssm 



~ /X {Hu) ~ 10^ GeV^ (A.8) 

Thus, Ojjt^jy^ generates a small Dirac mass term in a potentially viable range, while Oh^lNr 
generates a Dirac mass which is too big. 

It is also possible to consider operators which are quartic in the relevant superfields. As 
before, we restrict attention to operators which contain a factor of the form L^Nj^. There 
are precisely four possible quartic operators involving Hu, Hd, L^, Nj^ and X: 

^.^ XH UNj, ^ ^^^^ ^ 

^.^ X^H L^Nj, ^ ^^^^ ^ ^^^^^^ 
d^e ^^'f^'^ =^ qPQ (Nn) = +1, (A.ll) 



^X^HlUNi 



d'^ A =^ IPQ (Nr) = -7. (A.12) 

At 



In particular, only the second and third possibilities are compatible with an Eq GUT 
structure because the decomposition of the 27, 27 and 78 only contain U{1)pq of charges 
magnitude between zero and four. The estimated size of the Dirac mass in these two cases 
is: 

Auv ^uv 

^aqXH^^ ^ ^ 3. . ||k ^ 4 X 10-5 eV, (A.14) 

where in the ffist line we have used the fact that in F-theory GUTs, the Giudice-Masiero 
operator X^HuHd/Auy generates the /x-term in the effective theory. Both of these values 
are close to the required values for the neutrinos, although the ffist possibihty is somewhat 
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closer to the required mass scale necessary for matching to the observed mass splittings. 

Although it is in principle possible to consider operators with a larger number of fields, 
note that the largest field vev is set at the scale x ~ 10^^ GeV. As a consequence, each 
successive operator will be suppressed by a factor of roughly x/Mqut ~ 10~^, so that only 
quartic or lower operators are relevant for the present discussion. 



B Quartic Operator Dirac Scenario 

In this Appendix we discuss geometric configurations in F-theory GUTs which realize the 



quartic operators of lines (A. 10) and (A. 11) in Appendix A: 



As a piece of notation, we shall denote the triple intersection of three matter curves S^, 
and by S^S^Sc. Further, we shall often be interested in two configurations of triple 
intersections which share a common curve. For example, if the points joining HaHb^c 
and Sc-S/jS^; both lie on the curve Sc, we shall sometimes denote such a configuration as 
SaS^Sc ©Sc ^c'^d^e- We now turn to an analysis of various matter curve configurations 
which can generate the appropriate higher dimension operators. 

To see how the quartic operators are generated, consider a configuration where the X, 
Hu and H^, curves form a triple intersection such that X^H^Ha is gauge invariant such that 

also forms a triple intersection with the L and Nji curves. In terms of the same abstract 
ten-dimensional formulation provided earlier, the relevant interaction terms are given as: 



Leff D / d^eX^X + I d^eX^dAX + [ d^eUddAUd (B.2) 

JBi Jbz JB-i 

+ I d^ex'Hund + [ d^ew^cMR + h.c. (b.3) 

JBj, Jbs 

The first two F-terms originate from the covariant derivative on the appropriate curve. 

The second two F-terms originate from the triple overlap of matter curves. The X'^ and "H^ 
F-term equations of motion therefore contain the terms: 

^ Ud = Ha-^ {CUr) , (B.4) 



dm dA 

A-^X = X-^ {HJia) , (B.5) 
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so that: 

x = x- —nu— {i^J^r) + - (B.6) 

where we have dropped terms which will not figure in our discussion. Plugging this into 
the resulting D-term X'^X ^ we therefore obtain: 

L^ffDfdWH^UN'^.f Wx-lvl/^^-lvl/^^vl;^^. (B.7) 

In principle, this quartic operator can also be generated in configurations where L and 
do not even meet at a common point. In this case, we can consider a configuration where 
X, N and some additional curve meet at some point in the threefold base. Assuming that 
this curve also forms a triple intersection with Hu and L, it follows that by integrating out 
the massive modes S(BS'^ localized on the singlet curve, an analogous expression will again 
be generated. To be explicit, in this case, we consider the configuration of matter curves 
SxSatSs T^sT^lT^h^. We can write down the superpotential terms as before, solve the 
F-term equations of motion for X'^ and S'^, and plug the result back into the X'^X D-term. 
A similar analysis then yields the coupling: 

L^ff^ UwH^UNi- f ^^^^m^^m^mj^. (b.s) 

Finally, although we will not present explicit geometric configurations here, we note that 



the second quartic operator of line (B.l) can also be generated by integrating out massive 
modes. 



C Other Neutrino Scenarios 



In this section we collect some other possible neutrino scenarios which it would be interesting 
to develop further. Our aim here is not so much to provide an exhaustive list of alternative 
scenarios, but rather, to present some other potential avenues of investigation. To this 
end, we first discuss some additional Dirac mass scenarios where right-handed neutrinos 
localize in the bulk, and also discuss the numerology of instanton induced Dirac mass terms. 
After this, we briefiy mention another seesaw of potential interest based on massive string 
excitations. 



C.l Miscellaneous Dirac Scenarios 

In section [6] and Appendices A and B, we presented an analysis of Dirac mass terms where 
the right-handed neutrinos localize on curves normal to the GUT seven-brane. In that 
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context, higher dimension operators generated the necessary suppression in the mass scale 
below the scale of electroweak symmetry breaking. Here, we discuss scenarios where the 
right-handed neutrinos propagate in the bulk of the GUT seven-brane, and models where 
instanton effects can potentially generate a viable mass term. 



C.1.1 Nr From the Bulk 

So far in this paper we focused on the case where Nr lives on a matter curve. In this 
subsection we briefly note that in the Dirac scenario, it is also possible to consider models 



where Nr propagates in the bulk, This is especially natural in configurations where the 
bulk gauge group is of the form Eq, so that Nr descends from a spinor of SO (10) with PQ 
charge —3. On the other hand, as noted in f5,9|, such models typically contain extraneous 
zero mode states beyond those present in the MSSM. 

Putting aside this potential issue, we now consider a geometry where there is a local 
enhancement from Eq to Ej along curves, and Eq to Eg so that the corresponding modes 
trapped on the curves can form the 27^ interaction term. We consider configurations where 
X descends from the 27, while H^, L descend from the 27 and Nr from the 78. Note 
that the operator X'^ H^LNr is indeed invariant. Since is in the 27 we have the fusion 
rule [55j: 

78 x 27 = 27 + 351 + 1728, (C.l) 

the interaction term 27 x 78 x 27 x 27 contains a 27^ term, and thus a singlet as well. 

In addition to this zero mode content, we will also keep track of the 16 Kaluza-Klein 
mode excitations on the X-curve, which we denote by S'^. In this case, the relevant inter- 
action terms are: 

Le// D f d^eX^X + [ d^9S'dxS+ [ A^OX^dj^X 

Jbs Jbz Jbs 

+ [ d^esAfRX^+ [ d^es^-HuC + h.c. (c.2) 

Jb3 

where in the above, the first two F-terms originate from the associated kinetic terms on the 
X-curve. The third F-term originates from a coupling between a bulk gauge field and two 
chiral fields localized on the same curve (an S'SE coupling, in the terminology of [5]), and 
the last originates from the triple intersection of three matter curves. 

We now proceed to integrate out the relevant Kaluza-Klein modes for S and X. We 
solve the F-term equations of motion for X'^ and S'^, and plug the result back in the D-term 
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Although in different settings, there are neutrino models with right handed neutrino as bulk neutrinos. 



The main advantage is that the geometry gives the desired small mass. See 54 
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X'^X to obtain the operator: 

L^ff^ U'eX^H^UNi- f ^^^^^H.^^l^^i. (C.3) 

We caution however that for Nr to hve in the bulk, we must have an Eq bulk gauge group. 
In this case, it is not clear whether it is possible to obtain a low energy spectrum completely 
free of exotic fields. It would be interesting to study the consequences of such a scenario 
in greater detail, and in particular to establish whether potentially problematic exotics can 
indeed be removed from such models. 



C.1.2 Instanton Induced Dirac Masses 

Instanton generated Dirac mass terms have been considered in intersecting brane models. 



for example in 56 . In the context of F-theory GUTs, the characteristic size of instanton 



effects is determined by the requirement that the instanton induced Polonyi term: 

j d?eFxX = MIq j d2^g(4)X (C.4) 

responsible for supersymmetry breaking generates a value of Fx = MpQq(^4^ of order: 

^pq9(4) ~ Fx ~ 10^^ GeVl (C.5) 

Here, MpQ denotes the characteristic mass scale of the Peccei-Quinn seven-brane, which we 
shall take to be roughly the GUT scale, and g(4) denotes the suppression factor associated 
with a D3-instanton wrapping the same surface as the Peccei-Quinn seven-brane. Note that 
since X has —4 units of PQ charge, g(4) will have +4 units of PQ charge. Returning to the 
operator: 

'd'en^LNn, (C.6) 



when Hu, L and have respective PQ charges —2, +1, +1, this operator is invariant under 
U{l)pQ, and so will not be generated by instanton effects. On the other hand, when Np 
has PQ charge —3, so that it descends from the 78 of E^, the resulting operator HuLNp 



will have PQ charge —4 and so can in principle be generated by instanton effects The 
resulting coefficient is then given as: 

g(4) 10-'' - 10-''- (C.7) 



^''Although it is tempting to consider instanton effects which directly generate the operator (HuL)'^ / Any 
in a Majorana scenario, note that instantons wih generate such operators with unviably small coefficients. 
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As a consequence, the resulting Dirac mass term will be of order: 

^Sri - m ■ vu - 10-^ - 10-6 eV, (C.8) 

which is slightly too small. In principle, however, such effects could be present and may 
generate additional subleading corrections. It would be interesting to evaluate the expected 
flavor hierarchy derived from estimating the overlap of instanton zero modes. 

C.2 Symmetric Representation Seesaw 

As a final possibility, we consider another Majorana scenario derived from fields transform- 
ing in two index symmetric representations (the 15 or 15) of SU{5). In terms of represen- 
tations of the SU{2)l factor of the Standard Model gauge group, these fields transform in 
the triplet of SU{2). This scenario then realizes the triplet seesaw mechanism. 

We consider a configuration of matter curves where the Higgs up self-intersects, and 
the lepton doublet curve self-intersects, such that both self-intersections form a triple inter- 
section with a curve E15 where six- dimensional fields transforming in the 15©15 of SU{5) 
localize. In terms of the notation introduced in Appendix B, this can be described as 
the matter curve configuration S^/^Sj^^Sis ©215 ^w^l^l- Letting Ni^ © iVfg denote a 
vector-like pair of matter fields localized on E15, the superpotential will contain the terms: 

W D H^H^Nl^ + LLN^^ + MN^^N^^, (C.9) 

which would realize a variant of the Kaluza-Klein seesaw. Note that fields transforming in 
the two index anti-symmetric representation of SU{5) would not couple to Hu- 

The resulting light neutrino masses are either difficult to accomodate within a GUT 
framework, or tend to be too small. The essential problem stems from the fact that the 
15 of SU{5) is a two index symmetric representation of SU{5), and so as a massless six- 
dimensional field necessarily localizes on a curve where SU{5) enhances to USp{10). This 
can be arrived at by noting that the 10 of SU{5) locahzes on a curve of SO{10) en- 
hancement. Unfortunately, USp{10) does not embed in £"8; such a configuration is thus 
somewhat counter to the notion of E-type structures, which have figured prominently in F- 
theory GUTs. On the other hand, as is well known in the context of perturbative orientifold 
constructions, when the massless sector consists of fields transforming in the 10 of SU{5), 
the first excited string state will transform in the 15 of SU{5). Although the analogue of 
the perturbative string states are not known in the present context, it is likely that some 
massive modes localized on a curve where SU (5) enhances to SO{10) will indeed transform 
in the 15 of SU{5). Since these fields correspond to the analogue of massive string modes, 
they are quite heavy, and as such, will tend to have a seesaw suppression scale which is far 
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too high. We therefore conclude that resulting light neutrino mass scale again tends to be 
slightly too small. Even so, it could nevertheless be of potential interest to study such a 
scenario in more detail. 



D Haar Measure and Mixing Angles 

In this section we review the parameterization of the Haar measure for 3x3 unitary 
matrices in terms of Euler angles. The Haar measure must be used to generate random 
unitary matrices. Using the parameterization in terms of Euler angles, we extract the 
probability distributions for the three neutrino mixing angles for random unitary matrices. 

There are various algorithms to generate random unitary matrices using the Haar mea- 
sure. A particularly simple one, using parameterizations of unitary matrices in terms of 



Euler angles, is explained in 57 . Any 3x3 unitary matrix U can be written as 



U = e^"^(''=^)(^23,fe,r723)^^'''^(^i3,V'i3,0)^(^''H^i2,V^i2,r/i2) 
where the E^'^'^^ are 3x3 unitary matrices with entries 

hi 



(D.l) 



COS OijQ'-'^^i 

COS 6ije~''^'^ 
sin OijB^^'^ 



for fc, / = 1, 2, 3 and k, I ^ i,j 
for k = I = i 
for k = I = j 
for k = i and / = j 



(D.2) 



— sin 6ije for k = j and / 
The three angles and six phases take values in the intervals 

TT 



< 9ii < 



< tlJij,r]ij,a < 271. 



(D.3) 



To make contact with our parameterization of the neutrino mixing matrix given in (2.7), we 
note that the three angles % correspond to the three mixing angles. Out of the six phases, 
three are physically irrelevant, and the three other ones correspond to the CP violating 
phases S, ai and a2, where the latter two are only physical in the Majorana scenario. 

In this parameterization the Haar measure can be written down explicitly. Following 
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57 , it reads 



PuidU) = Cda II dr/,_i,, J] d^,,d (cos^^^'^*) %) , (D.4) 

l<i<3 l<i<i<3 

where C is some normalization constant. From this exphcit expression for the Haar measure 
we can generate random unitary matrices as follows. First, we generate random phases a, rjij 
and ipij uniformly distributed between and 27r. However, we must not generate uniformly 
distributed random angles 6ij. Rather, we first generate random parameters ^ij uniformly 
distributed between and 1. Then, the angles are given by 

6'i2 = arccos ^^]^2^ j , 6*23 = arccos ^^23^) 5 6'i3 = arccos ^.^Jg^ j . (D-5) 

In other words, the probability density functions for the theta angles that must be used to 
generate random unitary matrices are given by the functions: 

P(^^i2) = 2sin(^i2)cos(^i2), 

P(^^23) =2sin(^23)cOs(^23), 

P(ei3)=4sin(Mcos=^(^i3). (D.6) 
As a further check, we also generated 100,000 random unitary matrices using the numerical 



algorithm presented in 58 . We then extracted the mixing angles from these matrices, and 



indeed obtained the probability densities (D.6). 



^^The minor difference between our expression for the Haar measure and the one presented in 57 can 
be traced back to a different ordering in the product of the matrices iJ^*'-?' above. We use the ordering that 
makes contact with the standard parametrization of the neutrino mixing matrix (2.7 1. 
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